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Abstract 

We review vibrational and electronic properties of single and a few layer M0S2 relevant to 
understand their resonant and non-resonant Raman scattering results. In particular, the optical 
modes and low frequency shear and layer breathing modes show significant dependence on the 
number of M0S2 layers. Further, the electron doping of the M0S2 single layer achieved using 
top-gating in a field effect transistor renormalizes the two optical modes A] g and differently 
due to symmetry-dependent electron-phonon coupling. The issues related to carrier mobility, 
the Schottky barrier at the MoS2-metal contact pads and the modifications of the dielectric 
environment are addressed. The direct optical transitions for single layer-MoS2 involve two 
excitons at K-point in the Brillouin zone and their stability with temperature and pressure has 
been reviewed. Finally, the Fermi-level dependence of spectral shift for a quasiparticle, called 
trion, has been discussed. 
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1. INTRODUCTION 


Scientific challenges and potential applications of single and a few layer graphene in the 
last few years have given imputes to look for similar opportunities in other layered ma¬ 
terials like BN, M 0 S 2 , MoSe 2 and WS 2 . The field-effect transistors made of two-dimensional 
semi-conducting materials are of particular interest for technological applications. In past 
few-years, the enormously studied single-layer graphene has made its importance in sensors 
and high-frequency devices mn because of the high carrier mobility (~ 10 5 for suspended 
graphene [9]) and the response to the surrounding charges; although, in logical devices it 
is not very useful due to the vanishing electronic band gap at the K-points of the Brillouin 
zone (BZ). Recently, a single layer transition dichalcogenide material Molybdenum disulfide 
(M 0 S 2 ) held effect transistor has been demonstrated at room temperature to have a high 
on-off ratio ~ 10 8 [ID] , whereas ~ 10 11111 thick M 0 S 2 device gives a phonon limited-mobility 
value as high as ~ 700 cm 2 /V-sec [TT] . Production of single layer M 0 S 2 is not limited to the 
mechanical exfoliation method pioneered by Geirn et al for graphene, it is being produced by 
the chemical vapor deposition (CVD) [12], high temperature annealing method [13], solution 
method (HHl6j, intercalation method HM9] and so many other techniques [20H23] . CVD 
based M 0 S 2 is used with graphene to form heterojunction which produces photoresponsitiv- 
ity value ~ 10' A/W [T2]. MoS 2 nano-particles deposited on reduced graphene oxide act as 
an advanced catalyst for the hydrogen evolution reaction m The diamagnetic single layer 
M 0 S 2 shows a magnetic ordering at room temperature upon exposure of proton of 2 Mev 
energy which is attributed to the formation of various defects and edge states [25[ 126IJ . 

The photoluminescence (PL) spectra confirms the gradual transition from an indirect 
band gap of bulk-MoS 2 to the direct band gap of single la.yer-MoS 2 m- The presence of the 
spin-filtered valence band at the K-point valleys of BZ and the time reversal symmetry guar¬ 
antees the strong couplings of the spin and the valley index of the charge carriers and hence 
consequent optical selection rules [28] . It has been shown with the circularly polarized light 
that the observed chirality of the PL emission has almost the same degree of polarization as 
the incident one at low temperatures [29]. The voltage dependent transfer of the spectral 
weight ra of the exciton complexes ~ 1.9 eV in PL spectra is associated with the observed 
metal-insulator transition [31] in the two-dimensional M 0 S 2 at low temperature. The mea¬ 
sured ratio of the Coulomb potential energy to the kinetic energy is ~ 60 at a charge density 
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of ~ 1 x 10 11 cm -2 for the electron gas [30). Such high electron-electron interactions set the 
single layer M 0 S 2 as a platform to observe Wigner crystallization in the insulating regime, 
yet to be experimentally verified. Moreover, the existence of non-zero Berry curvature pre¬ 
dicts the Hall effect in presence of circularly polarized light without applying magnetic held 
[32j. The single layer M 0 S 2 shows an electric-held induced superconductivity in the metallic 
phase with transition temperature of T c ~ 9.4K [33]. This observed superconductivity is 
attributed to the electron-electron interactions and the weak electron-phonon interactions 
[31]. In this review, we will discuss various aspects of phonon-assisted phenomena and their 
probing by Raman spectroscopy as a function of temperature, pressure, carrier doping and 
the layer numbers; and the dependence of carrier’s mobility on various factors. Finally, the 
PL spectra at different gate voltages has been discussed leading to the discovery of negatively 
charged excitons or trions in single layer M 0 S 2 . 


2. CRYSTAL LATTICE AND ELECTRONIC BAND STRUCTURE OF MOS 2 : 
INDIRECT TO DIRECT BAND GAP 

Bulk Molybdenum disulphide (M 0 S 2 ) is a layered semiconducting material, each layer 
having hexagonal symmetry as graphene. Unlike graphene which consists of only carbon 
atoms, each M 0 S 2 layer unit is formed by stacking of three hexagonal planes of S-Mo-S 
atoms in order along the c-axis [35j as shown in Fig. [I] Within stable unit (S-Mo-S), atoms 
are bonded in a prismatic configuration through ionic-covalent interactions and each such 
stable units are connected via van der Waals interactions along the c-axis to form the bulk 
M 0 S 2 crystal. Two units, forming the unit cell of bulk M 0 S 2 , are stacked vertically in such 
a way that the Mo atom of one layer sits on top of the S atoms of other layer as shown 
in Fig. [I] There are two polytypes [36] of M 0 S 2 : (i) 2 H-M 0 S 2 [shown in Fig. [I] is more 
commonly available in nature as molybdenite salt and it has two layers per unit cell. Our 
discussions will be related to this type of structure; and (ii) 3R-MoS 2 has three layers per 
unit cell, which are stacked in rhombohedral symmetry. 

Like Si, bulk M 0 S 2 has an indirect band gap of ~ 1.2 eV |27[ 37, [38] with valence band 
maximum at T point and conduction band minimum in between the T — K direction as 
shown in Fig. [2j[a) . Now if we cleave the bulk M 0 S 2 one by one stable unit (S-Mo-S) and 
examine the electronic band structure, we will see that there is a transition from indirect 
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FIG. 1 . (color online)-The trigonal prismatic configurations of Mo atoms are shown for the struc¬ 
ture of 2H-M0S2 type. The (a,b) axis is in the x-y plane and the c-axis is along the z-direction. 
Along c-axis S-Mo-S units are stacked via van der Waals interactions. Taken from ref. j 35 l. 


band gap to a direct band gap ~ 1.9 eV at K point for single layer (S-Mo-S unit) M 0 S 2 
(SL-M 0 S 2 ). The calculated electronic band structure with the number of layers [3H] is shown 
in Fig. [3j The photoluminescence spectra [SUES] confirms the direct band gap of SL-M 0 S 2 , 
which will be discussed later. For the bulk or a few layers of M 0 S 2 , it has been shown m 
that the states near the indirect band gap consist of linear combinations of d orbitals on Mo 
atoms and p z orbitals on S atoms. The states near the conduction band at K point consist of 
only localized d orbitals on Mo atoms. Since Mo atoms are sandwiched between two planes 
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of S atoms in single S-Mo-S unit, the direct band gap has less dependence than the indirect 
band gap on number of layers (or interlayer coupling). The quantum confinement effect 
along the c-axis increases the indirect band gap, whereas direct band gap does not vary much. 



FIG. 2 . (color online)-Calculated electronic energy band-diagrams are shown for (a) bulk-MoS2, (b) 
4 -layer M0S2, (c) 2 -layer M0S2 and (d) single-layer MoS2- The solid arrow lines shown in Figs.(a), 
(b) and (c) indicate the indirect transitions from valence band (at T point) to the conduction band 
minimum (along T-K direction), whereas in (d) it indicates the direct electronic transition at the 
K-point. Taken from ref. [ 39 ]. 


3. PHONONS IN BULK AND SINGLE LAYER MOS 2 

Bulk M 0 S 2 belongs to the point group D 6 h (space group P3ml) pTTj . Each unit cell 
with lattice parameters a=3.12 A and c=12.3 A BH contains two molecular units shown 
in Fig. i m and hence, it has total 18 normal modes of vibrations corresponding to 
2A 2u +2E lu +2B2 ff +2E2 9 +A l9 -l-E lff +B lu -l-E2 U irreducible representations pf3l - H5] . The cal- 
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FIG. 3. (color online)-Calculated changes in the band gap for M 0 S 2 flakes as a function of (a) 
number (N) of layers and (b) the separation (w) between two-layers. Here A and B denote S-Mo-S 
unit. Taken from Ref. [3D]. 

culated values of the frequencies and optical activity of all the phonon modes are shown in 
Table 1 [0T] . K\ u and E] u phonon modes belong to acoustic branches. Since bulk MoS 2 has 
center of inversion symmetry, the infra-red (IR) active phonon modes (A§ u and Ef u ) and 
Raman active phonon modes are mutually exclusive. Excluding the inactive phonon modes 
(2B 29 , Bi u and E 2u ), the total four Raman-active optical phonon modes are A± g , Ei 9 and 
2E 2ff . The vibrational configurations for all the phonon modes are shown in Fig. [5|a) and 
the phonon dispersion curves for bulk MoS 2 are shown in Fig. [5](b) [12], 

Single layer MoS 2 has a point group symmetry of D 3 h (space group P6m2) [01]. Since 
SL-MoS 2 has no inversion symmetry, the labellings for the phonon modes gets changed. The 
total number of phonon branches is reduced to nine and the modes at T-point are given by 
2A2+2E , -l-A , 1 -l-E ,, [05]. Here Ag and E' are acoustic phonons; another E* mode is both IR 
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FIG. 4 . (color online)-Hexagonal lattice structure of 2H-M0S2 is shown in (a) for side and top 
views. The unit cell is marked by solid black lines, which contains one Mo atom (shown by red 
balls) and one S atom (shown by grey balls). The corresponding Brillouin zone is shown in (b), 
where red dots mark the high-symmetry points and dashed lines the corresponding directions. 
Taken from Ref. [ 12 ]. 


and Raman active. Second Ag mode is IR active. The total Raman active phonon modes 
are A], E* and E” corresponding to bulk Raman modes A ]g , E 2g and Ei g respectively. The 
eigenvectors of the phonon modes and the calculated dispersions are shown in Fig. [6](a) and 
Fig. |6](b), respectively [42j. Even number of S-Mo-S units belong to D^. (space group P3ml) 
point group symmetry having inversion center and odd numbers belong to D 3 h (space group 
P6rn2) point group symmetry without inversion center |41| . 


4. LAYER DEPENDENCE OF OPTICAL RAMAN MODES: OFF-RESONANCE 

With the ability of making single and bi-layer samples by mechanical exfoliation tech¬ 
nique, the layered materials have potential technological applications in nano devices and 
logic circuits. The non-destructive Raman spectroscopy has been widely used to measure 
the number of layers precisely for single, bi-layer and multi-layer graphene systems via the 
signatures of the 2D phonon band near ~ 2700 cm -1 (two-phonon double resonance feature) 
[48]. Similarly, different number of M 0 S 2 layers have been identihed through the mode fre¬ 
quencies. In this section, we will discuss the off-resonance Raman spectra (laser excitation 
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FIG. 5 . (color online)-(a) Different optical phonon modes at the T point with the corresponding 
symmetries and vibrational configurations for bulk-MoS2- The numbers are the corresponding 
vibrational frequencies in cm -1 . The phonon dispersion curves along different symmetry directions 
and the corresponding density of states are shown in (b) and (c), respectively. Taken from Ref. 
021 - 


energy 2.41 eV and 2.33 eV as compared to the direct band gap of 1.9 eV). 
















TABLE I. List of high-frequency optical phonon modes for bulk and SL-MoS2. 


SL-MoS 2 

Bulk 

Character 

Direction 

Atoms 


uj(cm 1 ) 

uj(cm x ) 

D 3 /,(S) 

De/i(E>) 



Involved 


calculated® 

exp. 

A" 

a 2u 

Acoustic 

Out of plane 

Mo+S 

0 

.0 (S) 0.0 

(B) 



B l 9 

Inactive 

Out of plane 

Mo+S 


55.7 

(B) 



E l 9 

Raman 

In plane 

Mo+S 


35.2 

(B) 

22 (Bi) 33 (B) c ’ d 

a; 

A 1 g 

Raman 

Out of plane 

S 

410.3 

(S) 412.0 

(B) 

402 (S) 408(B) b 


Blu 

Inactive 

Out of plane 

S 


407.8 

(B) 


a" 

a 2 

A2u 

Inactive (E 

11 c) Out of plane 

Mo+S 

476.0 

(S) 469.4 

(B) 



b 1 9 

Inactive 

Out of plane 

Mo+S 


473.2 

(B) 


e' 2 

E 2 9 

Raman 

In plane 

Mo+S 

391.7 

(S) 387.8 

(B) 

382 (S) 380(B) b 


Eiu 

Infrared (E 

T c) In plane 

Mo+S 


391.2 

(B) 


e 2 

E lg 

Raman 

In plane 

S 

289.2 

(S) 288.7 

(B) 



E^u 

Inactive 

In plane 

S 


287.1 

(B) 



a Ref. [fij. b Ref. [3B] , c Ref. [50] . d Ref. |54j . Here ’Bi’ indicates Bi-layer 


4.1 . High frequency optical modes (around 400 cm x ): opposite trends of A\ g 
and Bl g modes as a function of number of layers 

The single and few layers of MoS 2 deposited on Si0 2 /Si wafer are shown in Fig. [7](a) 
taken by optical microscope and the thickness of the corresponding flakes measured by 
atomic force microscopy (AFM) are shown in Fig. [7](b). The measured thickness of a single 
layer MoS 2 is about 0.6-0.7 nm. The Raman spectra for bulk and few layers MoS 2 are shown 
in Fig. ||(a) and the observed frequencies of the two modes are shown in Fig. [8](b) [4?]. The 
Ai 9 Raman mode hardens with increasing the number of layers by ~ 6 cm -1 , whereas 
softens by ~ 2 cm” 1 . Above four layers, both modes show saturation values with the Bulk 
values ~ 408 and 382 cm” 1 , respectively. The opposite trend of the two Raman modes 
leads to increase in frequency difference [shown in Fig. [8|c)] which can be used to probe 
the number of layers of MoS 2 system. The hardening behavior of out-of-plane mode (Ai g ) 
can be understood in terms of increasing the restoring force constant perpendicular to the 
basal plane through van der Waals interactions as we increase the number of layers. In 
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FIG. 6. (color online)-(a) Different optical phonon modes at the T point with the corresponding 
symmetries and vibrational configurations for single layer-MoS 2 - The numbers are the corre¬ 
sponding vibrational frequencies in cm -1 , (b) Phonon dispersion curves along different symmetry 
directions and (c) the density of states. Taken from Ref. [42] , 


comparison, the in-plane vibrational mode (E^) shows anomaly. The difference between 
these two vibrational modes is that involves vibration of Mo and S atoms in the basal 
plane; whereas in A Ug vibrational configuration, Mo atoms remain fixed and only S atoms 
vibrate along c-axis. It has been shown 05 that the long range coulomb interaction part 
of the self interaction remains negligible for S atoms, whereas for Mo atoms it decreases 
considerably as we increase the number of stacking layers. Because of increased screening 
and consequent increase in dielectric tensor with the number of layers, the in-plane mode 
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FIG. 7. (color online)-(a) Optical image of the M 0 S 2 flake deposited on Si02/Si wafer. The dashed 
box indicates the area of 8 x 8 /i m 2 for which the AFM image has been taken and shown in (b). 
The different layers and the corresponding height profiles are shown in AFM image. Taken from 
Ref. [35] , 


El g softens. The qualitative agreement between the theory [JIJ and the experiment [46] for 
the observed two Raman modes is shown in Fig. [9} 


The integrated intensity variation of the Ai ff and modes with the number of layers 
is shown in Fig. [Tdj / a). The Raman intensity increases linearly upto four-layers and then 
decreases to the bulk value. The increase in intensity for both the Raman modes can be 
attributed to the optical field enhancement due to multiple reflections of Raman field within 
the sample layers, and consequent interference effects between these reflected lights and the 
those coming from the thin layer of SiC >2 (300nm) deposited on Si substrate. The optical field 
enhancement and interference phenomena also leads to increase in intensity of Raman modes 
for graphene systems and shows a maximum upto 9 layers [3Hj. It has been observed [H] 
that SL-M 0 S 2 deposited on quartz substrate has Raman intensity of ~ 20 % of that of the 
SL-M 0 S 2 deposited on SiC^/Si substrate. The integrated intensity ratio of the two modes, 
I(E 2 S )/I(Ai 9 ), is also shown in Fig. 10 a). Since the interference enhancement phenomena 


affects both the modes equally, the abrupt increase of the integrated intensity ratio from 
single to bi-layer is not clear and, it may be due to the differences of underlying crystal 


11 




(a) 



Raman shift (cm’ 1 ) 


(b) 



(c) 



FIG. 8. (color online)-(a) Raman spectra of a few layers and bulk M 0 S 2 . The laser excitation 
energy used is 2.41 eV (514 nm) 0J]- Solid lines represent the Lorentzian fits to the experimental 
data (circles), (b) Raman shifts of these two modes and (c) the frequency difference {uA la -w E i ) 

9 2 g 

as a function of thickness of the sample. Taken from Ref. 03 - 


symmetries between them. The full width at half maximum (FWHM) of A\ g and modes 
are shown in Fig. 10 b) 07]. The FWHM of A lg Raman mode increases from 5 to 7 cm -1 
as layer number decreases from 7 to 2 layers. The E\ mode shows no thickness dependent 
variations (~ 1 cm -1 ) of the FWHM. The FWHM dependence of both the modes A\ g and 
E l g needs further understanding. Lee at el 06] showed that the broadening of the Ai g 
phonon is also observed for the suspended M 0 S 2 layers and hence cannot be attributed to 
the interaction with the substrate. 
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FIG. 9. (color online)- The phonon frequency differences between n-layers and bulk for A\ g (top 
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FIG. 10. (color online)-(a) The integrated intensity of A\ g and modes (corresponding to left 
vertical scale) and their ratio (corresponding to right vertical scale). Taken from Ref. [56]. (b) 
The line width (FWHM) of these two modes. Taken from Ref. [57] , 


4.2 . Low frequency optical modes (below 50 cm 1 ): shear (E \ g ) and compressional 
modes 


In the E\ g optical phonon, which is a in-plane rigid layer vibration (Fig. Ill), the atoms in 
each layer vibrates in phase such that each layer rigidly oscillates against the adjacent layers 
and the only restoring force is the weak van der Waals interactions acting between the layers. 
This mode, also known as a shear mode, exists in all layered compound materials with low 
frequency (~ 30-40 cm -1 ). Raman spectra [50] of a few layers and bulk M 0 S 2 for the low 
frequency optical phonons are shown in Fig. |TT|[b). The E\ g mode shows stiffening behavior 


with the thickness of the sample from ~ 22 cm -1 (2 layers) to ~ 33 cm -1 (bulk) [Fig. 11 3 ]. 
Another low frequency mode with broad feature attributed to the compressional mode (C- 
mode) shows softening behavior from ~ 42 cm -1 (2 layers) to ~ 15 cm -1 (7-layers). The 
experimentally observed and theoretically calculated [50j frequencies for both the modes 
are shown in Fig. 11(c). For the shear E\ g mode, the fitted frequency (Fig. 11 3 ) follows the 


expression (obtained from the linear-chain model) u> = + cos(j^) [5Tj. Here cu 

denotes the peak position of the shear mode; N, the number of layers; c, the speed of light 
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in cm/s; a, the interlayer force constant per unit area and / x , the unit layer mass per unit 
area. For multi-layer graphene px[], the low frequency shear modes also follow the same 
trends i.e. with the increasing layer numbers, the restoring force constant acting between 
the layers increases. 
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FIG. 11. (color online)-(a) The vibrational configurations of the E \ g and the compressional modes, 
(b) Raman spectra of two modes for different layers and bulk M 0 S 2 deposited on Si02/Si substrate. 
The dashed line is guide to the eye. (c) The peak positions as a function of thickness of the samples. 
The interlayer force constant per unit area (a) and the unit layer mass per unit area (p) used for 
the least square fitting of E 2 g mode are 2.9 x 10 19 N m -3 and 30.3 Kg A~ 2 , respectively. Taken 
from Ref. ESDI 
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The thickness dependence of the C-mode shows ^ dependence instead of The A 
behavior also observed for the longitudinal acoustic phonons on Na films [deposited on Cu 
(001)] [22], is due to open standing waves in multilayer (or thin films) systems [22]- A 
strong coupling between the substrate and the atomic layer of the sample does not allow the 
atoms close to the substrate vibrate in the C-mode and creates a node there, while there is 
no restrictions for the atoms in the top layer; thus creating open standing waves (kind of 
organ-pipe modes) and hence A dependence. 


4.3 . Thickness monitoring of M 0 S 2 flakes via Raman imaging: two pairs of Raman 
modes 


From the above two subsections, we have seen that the frequencies of the four Raman 
active modes (2 -pair) change differently by sensing only the number of M 0 S 2 layers. The 
2-pairs of optical phonons are {A ls and E\ g } and {E \ and C-mode}. Each member of the 
first pair shows opposite trend with the result that the difference u>A lg - oj e * shows a sizable 
change, which helps in Raman imaging. Since, the C-mode shows softening behavior and it 
has low frequency (15 cm -1 for 7-layers) than the E \ g mode (32 cm -1 for 7-layers), the shear 
E 2 2g mode is the feasible one from experimental point of view. Moreover, the shear mode 
shows a ~ 50 % changes of peak positions going form bi-layer (22 cm -1 ) to bulk value (33 


cm 


— 1 > 


Plechinger et al recorded (shown in Fig. fL2J) Raman images of a few layers of 


M 0 S 2 by locking the frequency difference of two high frequency Raman modes (u)A la - Wgi ) 

9 2 g ' 


and a low frequency shear mode (E \ g ). The Raman imaging of M 0 S 2 flakes (rectangular 


area marked in Fig. 12 1 ) is shown in Fig. 12 b) for the frequency difference of ua 1 -oj e 1 • 

_ _ ^ 2 g 


The single layer M 0 S 2 for which, coa 19 ~ =18 cm 1 (from Fig. 8p), is shown by marked 


area with false color in Fig. 12 b). The signature of E\ g mode on single layer is zero and is 
shown in Fig. [l2|(c) by a closed marked region. The bi-layer (2L), tri-layer (3L), hexa-layer 
(6L) and 10-layers (10L) and their boundaries are clearly seen by tracking the shear mode 
(in Fig. Il2fc), whereas the counterpart for the frequency difference is not clear as shown in 


Fig. 12 » 
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FIG. 12. (color online)-(a) The optical image of different layers of MoS 2 - The boxed area used 
for Raman imaging, (b) The Raman image taken by locking the L0diff= uja 19 -^e] • Single layer is 
marked by a boxed area, (c) Raman image recorded by locking the uj s h ear =oj E 2 . Here the color 
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scales used for Udi // and uishear are arbitrary. Taken from Ref. [51] • 

4.4 . Shear and layer breathing modes of N-layers of M 0 S 2 { N = 1 —> 19 } 

For the rigid layer vibrations of N-layer M 0 S 2 , there are N-l shear modes parallel to 
the basal plane (perpendicular to the c-axis) and N-l layer breathing inodes (LBMs) along 
the c-axis. For even N, LMBs are not Raman active and A shear modes are Raman active 
(doubly degenerate). For odd, LBMs and N-l shear modes are Raman active jJ5]. 
Fig. |l3|(a) shows the newly observed low frequency shear and LBMs Raman modes for 5L 
and 6L with different polarizations p[5]. As stated earlier, the Raman active modes A) 
and E* for odd-layers (ONL-M 0 S 2 ) are the counterpart of the A\ g and E| for even-layers 
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(ENL-M 0 S 2 ), respectively. The eigenmodes of A\ and are shown in Fig. 13 b). 
mode is Raman inactive and the reason of observation is discussed later. All the 
low frequency modes are listed in Table If. 
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FIG. 13. (color online)-(a) Polarized Raman spectra for 5L and 6L of M 0 S 2 . YX represents the 
direction of incident laser polarization (Y) and the direction of analyzer’s polarization (X); both 
lie in the basal plane, (b) The vibrational configurations of the atoms in and B \ g modes. Taken 
from Ref. [ l5| . 


To understand the low frequency rigid layer modes, the calculated 0 a frequencies with 
the diatomic chain model (DCM) and the experimentally observed ones for the ENL-M 0 S 2 


and ONL-M 0 S 2 are shown in Fig. 14, separately. For the ONL-M 0 S 2 , the shear mode E 
(28 cm -1 for 3L) stiffens and the one LBM A', (29 cm _1 for 3L) softens with increasing 


the thickness of the sample from 3L to 19L (see Figs. 14 b and d). According to DCM 

calculations, a new mode appears as N reaches a value of 4N+3, where N=l,2,3,. And 

that new mode splits into two branches such that one branch softens and other stiffens until 
N reaches a value of 4N+3 and so on. The similar softening and stiffening behavior of a 
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TABLE II. List of low-frequency optical phonon modes for 5-layer and 6-layers of M 0 S 2 [35]. 


5-Layers (5L) 

6-Layers (6L) 

Shear modes 

Exp. 

LBMs 

Exp. 

Shear modes 

Exp. 

LBMs 

Exp. 

2E' (R, IR) 

19 & 33 a 

2a; (r) 

17 a & 47 b 

(R) 

23“ & 32 b 

3B| S C 

15 & 41“ 

{XX, YX} 

(in cm" 1 ) 

{XA} 

(in cm"" 1 ) 

{XX, YX} 

(in cm -1 ) 

{XX} 

(in cm"" 1 ) 

2E" (R) 

- 

2A" (IR) 

- 

2Ei u (IR) 

- 

2AL(IR) 

- 

{XZ, YZ} 

- 

- 

- 

- 

- 

- 

- 


a Ref. [45] . b Ref. [50]. c This mode is a silent mode; it is neither Raman nor infrared (IR) active. 


shear mode (23 cm -1 for 2L) and a LBM mode B 2g (40 cm -1 for 2L) for ENL-MoS 2 are 
shown in Figs. 14 (a) and (b), respectively. For ENL-M 0 S 2 , the appearance and splitting 
of a new mode into two branches follows according to 4N+2. The DCM calculated value 
matches well with the experimentally observed Raman inactive B 2g mode and this is the 
reason, along with the polarization behavior that Zhang et al [45] named it as a B 2g mode. 
Although, the B 2g for ENL-M 0 S 2 follows exactly the same polarization behavior as A] for 
ONL-M 0 S 2 (see Fig. 13 1 ), the deeper understanding of its appearance, polarization behavior 
and matching with the DCM calculated value needs further work. The most striking features 


from Figs. [14] is that the experimentally observed shear modes are coming from the upper 
branch of cone like curve with stiffening character and LBMs are from the lower branch with 
softening feature for both the ENL-M 0 S 2 and ONL-M 0 S 2 . The ratio of the shear E| mode 
for bulk to 2L is given by u>Buik/^ 2 L=^- 7/22.6=1.447 [45], which is very close to the value 
of graphene (u Bu ik/u 2 L = V%) [51]. 

In general, in-order to know the exact trend of softening and stiffening features of low 
frequency rigid layer modes of layered materials, a reduced monatomic chain model (MCM) 
has been considered [15]. The DCM calculated values showed that the relative displacements 
between Mo atoms and S atoms are so less (~ 0.6% for 2L) [45] and it decreases with 
increasing N. Therefore, for rigid layer modes, the relative displacements can be ignored and 
hence, in MCM model the additive mass ( M a dd = Mmo + 2 Ms, instead of a reduced mass 
4 = + is taken. All the observed shear modes for ENL-M 0 S 2 (E| ff ) anc i ONL-M 0 S 2 

(E') from Figs. 14(a) and (b) are plotted together in Fig. 15(a) with the DCM and MCM 
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FIG. 14. (color online)-Calculated (red cross) according to the diatomic chain model (DCM) and 
the experimentally observed frequencies or wavenumbers (blue circle) (a) with the even (ENL) and 

(b) odd number (ONL) of layers for the shear modes, and for the layer breathing modes (LBM) in 

(c) and (d). Here the connected solid lines are guide to the eye. Taken from Ref. |45| . 


calculated values m ; and Fig. |l5[b) for all the LBMs. The interesting feature from Fig. h5 


is that all the experimentally observed low frequency modes originate from ENL-M 0 S 2 he. 
2L, 4L, 6L and 10L. The black solid lines and the gray-dashed lines are the fitted curves 
which, follows ui(N) = u(2N 0 )^/l ± cos(N 0 n/N) (+ and - sign for upper and lower branch, 
respectively ) such that N > 2N 0 and iVo=l,2,3,..etc. Since, u>(2N 0 ) is almost equal to uj( 2) 
for all rigid layer modes, the fitted frequency obeys uj(N) = u( 2)y/l ± cos(N 0 n/N ) formula 
according to MCM. Here, the model (MCM) did not take into account the interactions 
between one layer of M 0 S 2 and the supported substrate and explains well the experimentally 
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FIG. 15. (color online)-Calculated (red cross) according to the diatomic chain model (DCM) and 
the experimentally observed frequencies or wavenumbers (blue circle) of all the layers together 
for the shear modes in (a) and the LBM modes in (b). Size of the circle represents the observed 
Raman intensity. The black solid lines (in a and b) are fitted curves according to the monatomic 
chain model (MCM) for those modes generating from 2L (IVo=l), 6L (Nq=3) and 10L (IVo=5) 
branches. Another set of modes (shear and LBMs) according to the MCM model are shown by the 
gray dashed lines in (a) and (b). Taken from Ref. [15.]. 


observed frequencies; whereas in contrast, Zeng et al |50j included the substrate effect and 
showed the 1/N behavior of one LBM as shown in Fig.[IT|c). It was shown that for suspended 
multilayer graphene, the scaling of the shear mode with thickness follows the trend predicted 
by MCM calculations [51 j . Recently, suspended 2L-MoS 2 shows no substrate effect on the 
rigid layer phonon frequencies and their FWHMs |55| . It would, therefore, be reasonable to 
assume that the interaction between substrate and MoS 2 is not responsible for the observed 
1/N scaling with the thickness J35J. 


Recently, two new low frequency rigid layer modes are observed |55j (Fig. 16 1 ): one is 
LBM (B2) and another is a shear mode (S2). Both modes appear from 4L-MoS 2 such that 
B2 mode softens from 54 cm -1 (4L) to 35 cm -1 (7L) and S2 mode stiffens from 13 cm -1 
(4L) to 30 cm -1 (12L). The frequency variation with the thickness of the sample is shown 


in Fig. 16 b). The fitted solid curves are according to linear chain model calculations 
S2 mode shows softening behavior because, the out of phase displacement between layers 
increases with N [55]. B2 mode stiffens as in-phase motion increases with the thickness N. 
The appearance of B2 LBM mode for even layers of MoS 2 (4L) needs further studies as, 
according to Zhang et al [45] for even layers, there should be no Raman active LBM. 
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FIG. 16. (color online)-(a) Raman spectra of few layers and bulk M 0 S 2 for the low frequency 
shear (S) and layer breathing modes (B). Here 1TL means one S-Mo-S unit i.e. one layer (1L). 
The dashed lines are guide to the eye. (b) The variations of the measured (solid dots) and the 
calculated frequencies (open squares) from first-principles. Solid lines represent the fitting of the 
data according to the linear chain model. Taken from Ref. |55j . 


5. LAYER DEPENDENCE OF OPTICAL RAMAN MODES: RESONANCE 


Resonance Raman scattering (RRS) occurs for a system, when the incident laser excita¬ 
tion energy is close to the electronic band gap. I 11 this section, we will discuss the RRS with 
the laser energy 1.9 eV (633 nm), which matches with the direct band gap at the K-point of 


the BZ. Fig. 17 shows the resonant Raman spectra m of different layers (N=1L, 2L ,4L and 
7L) and bulk M 0 S 2 . Since, the RRS enhances the possibility of higher order Raman scatter¬ 
ing (multiphonon process) compared to the first order process in normal Raman scattering, 


Fig. 17 shows many Raman peaks along with the off-resonant Raman modes. Figs. [l8[a), 
(b) and (c) show the variations of frequencies, FWHMs and integrated intensity ratio of two 
Raman modes Ai g and E^, respectively for the RRS [37]; which show the same trend as 


for the off-resonant case (see Figs. [8p, and 10). The new Raman bands, which were not 
observed with the 514 nm and 533 nm laser lines, are understood as follows. 


22 








































FIG. 17. (color online)-Resonant Raman spectra of few layers and bulk M 0 S 2 using 633 nm laser 
line. The strong mode coming from the Si02/Si substrate (marked by *) is absent for thick samples 
(bulk). Here N represents the number of layers. Taken from Ref. [47] . 


5.1 . New modes in low frequency region of 160 - 230 cm 1 

(i) A Raman band ~ 179 cm -1 is observed and this mode is assigned as a difference 
combination mode of A lp (M)-LA(M) [551 [57J. The frequency of this mode does not show 
layer dependent feature. The three Raman active (Ai s , E 1 g and E^) phonons at T point of 
BZ are almost dispersionless along the T-M direction; the frequency of rigid layer mode E| 
increases with q (phonon wave-vector) and reaches to a longitudinal acoustic (LA) mode at 
M (233 cm -1 ) [5?]. At M point of BZ, the dispersion of LA mode shows an inverse parabolic 
nature, i.e. on either side of M-point, the frequency decreases. The asymmetric nature of 
the peak on the higher side is due to the involvement of the LA(M) mode as a difference 
combination band (571- (ii) Another band at ~ 233 cm 1 is observed for 2L, 4L and 7L 
and is assigned [58J as LA(M). The activation of this first order mode implies the relaxation 
of the q ~ 0 selection rule. The appearance of the zone boundary phonon (LA) for M 0 S 2 
nanoparticle is attributed to the defect induced Raman scattering [59]. In the present case, 
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FIG. 18. (a) Raman shifts and (b) the line widths (FWHM) of the A.\ g and the modes as a 
function of the number of layers, (c) The integrated intensity ratio I(Ai 9 )/I(E 2 S ). Taken from Ref. 

m- 

the observed LA(M) mode for 2L, 4L and 7L MoS 2 could be due to the structural defects 

m- 


5.2 . Frequency region of 360 - 430 cm 


-l 


There are three new Raman active bands in the frequency range of 360 - 430 cm 1 . 

377 cm -1 (shown in Fig. 17) becomes sharper with increasing 


(i) a band named as c 
the layer from 1L to 7L. Sekine et al [60J assigned this mode as Ef u mode. Since, the 
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FIG. 19. (a) Raman shifts of the 2LA and the A. 2 U modes as a function of the number of layers. 
(b)-(d) The changes of the integrated intensity ratios for different modes. Taken from Ref. [ 17] . 


experimentally observed mode has higher frequency than the mode [23* 61] [62] and 
these are Davydov pair [SU], the observed c mode cannot be assigned as a Ef u mode [42]. The 
c mode does not show thickness dependent frequency variations. The integrated intensity 


ratio I(c)/I(E 2 S ) decreases with the number of layers as shown in Fig. 19 b). (ii) One peak 
~ 409 cm -1 evolves with number of layers and is assigned as Bi u mode [EH [62]. The B lu 
and Ai 5 mode are the small frequency splittings of Davydov pair [[60] and the appearance of 
the B lu mode is due to the resonance effect. Fig. [l9tc) shows that the integrated intensity 
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ratio I(Bi„)/I(Ai s ) increases with the layers of the sample, (iii) The prominent feature at 
~ 420 cm -1 , marked as V in Fig. 17, is observed under resonant conditions. This mode 
is assigned as a two-phonon process and the appearance of this mode is related to the A 
exciton (~ 1.9 eV) at the K-point [581460] . The dispersion curve for polariton-exciton is 


shown in Fig. [20} Here, the intermediate polariton state involves the two phonons under 
resonant condition. First, the photon-like state from the high energy upper branch (inner 
branch) scatters to exciton-polariton state (outer branch) by scattering a longitudinal optical 
phonon (this dispersive quasi-acoustic phonon is the silent B\ g mode at the Y point) and 
this interaction is due to the deformation potentials; second, the polariton state decays to 
photon-like state by scattering the transverse optical phonon (with K 2 lu symmetry) [60] . This 
mode is a combination band, cob = ojqa + ^>to [121 • 


5.3 . Frequency region of 440 - 490 cm 1 

(i) The mode ~ 450 cm -1 (for bulk) is attributed to second order Raman (2LA(M)) 
[HUES ED], and (ii) the mode ~ 466 cm 1 (for bulk) is assigned as A 2u [59[E3J. The peak 
positions of the two modes (2LA and A 2 <J increases with the number of layers (N) as shown 
in Fig. [I9|[ a). Fig. |l9|(c) shows the decrease of the integrated intensity ratio l(2LA)/I(A 2u ) 
with the thickness mm- 


5.4 . Frequency region of 510 - 630 cm 


-l 


For the bulk MoS 2 , (i) the observed mode ~ 526 cm -1 is assigned as Ei 3 (M) +2 LA(M), 
(ii) ~ 571 cm -1 mode as 2E l9 (r), (iii) ~ 599 cm -1 as E\ g (M) + LA(M) and (iv) ~ 642 cm -1 
as A i g (M) + LA(M) [59], In all the multiphonon processes involving LA(M) mode, the sum 
combination band has asymmetric tail at lower side 152. All the new high frequency bands 


~ 554 cnr 1 (marked as Q in Fig. 17), ~ 588 cm -1 (marked as P) and ~ 628 cnr 1 (marked 
as S) are observed under the resonance condition for 2L, 4L and 7L MoS 2 [47] and were not 
observed for bulk MoS 2 and nanoparticles of MoS 2 [59] . The peak positions of these modes 


do not show variations with N. Fig. 19 [d) shows the variations of the relative integrated 


intensity of the bands (Q, P and S) with respect to the neighboring band. 

Fig. 21 shows another new broad feature ~ 38 cnr 1 (5 meV) [50] observed under the 
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FIG. 20. (color online)-Schematic diagram for the dispersion curve of a polariton and the two 
phonon Raman scattering involving a longitudinal quasi-acoustic (Q.A.) mode and the Ef u mode. 
The dashed and the solid lines marked by arrows represent the a and /3 processes, respectively. 
Here E^ and E^ represent the energy of Is level for the A and B excitons, respectively. Taken 
from Ref. [60]. 


resonant condition for the cross polarization geometry. The peak position and the shape of 
the mode do not change with the layer numbers. Since, the binding energy of the excitons 
(A and B with energies 42 and 134 meV, respectively [64]) are much greater than the 
observed peak with 5 meV, the appearance of this mode is not responsible for the exciton 
mechanisms. Defects and impurity effects are also excluded due to having energy comparable 
to the excitons in optical measurements [50] . The oscillating collective modes (plasmons, 
charge density waves) are also excluded due to the reluctant behavior of the position and 
shape of the peak; as the long-range Coulomb interaction would be screened with the layer 
thickness [50]. The spin-orbit coupling is responsible for the splitting near the top of the 
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FIG. 21. (color online)-Resonant Raman spectra using 633 nm laser line for samples of different 
thickness. The new broad feature ~ 38 end 1 appears only under resonant condition. Taken from 
Ref. [50] . 

valence bands at the K-point of BZ, which leads to the origin of A and B excitons. Zeng 
et al [50] argued that the same spin-orbit coupling mechanisms also splits the conduction 
band at the K-point of ~ 4 meV, which leads to the observed broad feature (~ 38 cm -1 ) 
associated with the resonant Raman scattering. The understanding of the appearance of all 
the new Raman modes under resonant condition requires further theoretical work. 


6. TRANSPORT CHARACTERISTICS AS A FUNCTION OF TEMPERATURE 
AND CARRIER CONCENTRATION: VARIOUS ASPECTS TOWARDS IMPROVE¬ 
MENT OF MOBILITY 

In technological applications of nanodevices, the two-dimensional (2D) layered-materials 
have a lot of potential due to easy fabrication. The highest mobility observed for the 
suspended graphene device is larger than 10 5 cm 2 / V-sec [9j. Due to the lack of an electronic 
band gap in single layer graphene, the 2D semiconducting materials like M 0 S 2 (on-off ratio is 
10 8 [Ij]]) family compounds has potential applications in digital electronics with low power 
dissipation. In order to replace silicon based logic-devices, a high current on-off ratio ~ 




10 4 to 10' [9] and the electronic band gap larger than 400 meV [65] are reasonable. The 
various scattering mechanisms (intrinsic and extrinsic effects), which are responsible for the 
momentum relaxations and thereby decreasing the mobility of charge carriers in M 0 S 2 , are 
discussed in the following sections. 


6.1 . Temperature dependence of carrier mobility in bulk-MoS 2 

Bulk M 0 S 2 is a layered material in which each individual 2 D layers are stacked along the 
z-direction via van der Waals interactions. To capture this anisotropic feature, a potential 
model was invoked [66] such that, total potential V(r)=V(x,y)+V(z); where V(x,y) and V(z) 
represent the potential in basal plane and perpendicular to the basal plane, respectively. 
According to this model, the excess carriers feel a series of parallel potential wells along the 
perpendicular direction and between each of them, there exist local energy levels for the 
carriers which depend upon the local width of the wells. This implies that the localization 
energy within each layer depends on the layer thickness. The variation of the layer thickness 
associated with thermal phonons may give a friction to the carriers. The energy of the excess 
carriers (due to decoupled potential) is given by E(p) = * y + 2 I z (cos^^- — 1) [67]; where 
rn x is the effective mass in basal plane, d z the inter-layer spacing, I z small overlap energy 
between adjacent layers and the second term of E(p) corresponds to the tight binding model. 
The corresponding density of states (DOS) is given by D(E ) = cog -1 (l — E/2I Z ) for 
0 < E < AI Z and D(E) = (pff for E > AI Z (67j; where N z denotes the number of 
layers/lcngth along the z-direction. The DOS is almost constant except for the low energy 
range, which is the characteristic feature for 2 -dimensional electron gas. 

In most isotropic semiconductors (non-polar), the high temperature mobility is limited 
by the quasi-elastic scattering with acoustic phonons via deformation potentials [68]. The 
acoustic phonon mode corresponds to local dilatations of the lattice. Since, between two 
adjacent layers there is a weak van der Waals force, the change in deformation potentials 
(CDP) due to the acoustic modes polarized in the basal plane will be more than that by the 
vertically polarized ones. Therefore, in layered semiconductors, the carriers will be scattered 
significantly by the horizontally polarized acoustic phonons. The calculated mobility due to 
the deformation potential scattering by acoustic modes is ^ T~ iac ° usttc [67]; where 

7 'acoustic _ ^ w hi c p i s quite different than the well known value of 7 acoustlc = 3/2 for isotropic 
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case because of peculiar behavior of DOS. 



Temperature [K] 


FIG. 22. Hall mobility versus temperature for the n-type bulk-MoS 2 - The solid line represents 
the fitting to the observed data. At low temperature, the increment of mobility is due to extrinsic 
properties. Taken from Refs. [66l [67] . 

In contrast to the acoustic mode, the optical inodes corresponds to the distortion of 
the unit cell. In isotropic materials (non-polar), the CDP due to optical modes is lit¬ 
tle and can be neglected. But in layered materials, the deformation potential is highly 
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anisotropic [67]. Those optical modes which can modulate the thickness of the individual 
layer will interact with the charge carriers as discussed above (because of series of parallel 
potentials along z). These thickness modulating optical modes which do not produce first 
order dipoles irrespective of having effective charges on each atom are known as homopolar 
modes. The mobility due to the vertically polarized (along z) homopolar optical modes is 
l-hmik" — Ho(T/T 0 )~' yhp ° Iar |[67j; where 7 hpolar > 1 and fi 0 = fj, bu i k at T=T 0 . In case of polar 
materials, where long range Coulomb interaction associated with polar optical modes cou¬ 
ples the charge carriers (via Frohlich interaction), it has been shown that r y polar = ryhpoiar _ ^ 


The experimentally observed Hall mobility of bulk M 0 S 2 is shown in Fig. 22 The 
fitted solid line corresponds to ii e bulk = 100(T/T 0 ) -2 ' 6 cm 2 /V sec, where T 0 =300K. To get a 
comparison of different scattering mechanisms of charge carriers by thermal phonons, the 


exponent ( 7 ) for the mobility has been plotted with the phonon energy as shown in Fig. 23 
Since, for bulk M 0 S 2 the only one homopolar mode ra corresponds to 60 meV (A lff ~ 408 


cm 1 for bulk), it is clear from Fig. 23 that the vertically polarized non-polar (homopolar 
mode) mode is responsible for the high-T limited mobility. 


6.2 . Temperature dependence of carrier mobility in single layer M 0 S 2 

From theoretical point of view: The electrical properties of single layer M 0 S 2 tran¬ 
sistor have been probed extensively in past few years. The different phonon scattering 
mechanisms with the charge carriers discussed in the previous section is calculated [US] and 
plotted at T=300K as a function of energy for nondegenerate carriers, as shown in Fig.[24|(a). 
The Brillouin zone for SL- M 0 S 2 is hexagonal and similar to graphene, two valleys exist at 
K and K' point. Kaasbjerg et al [69J calculated both the intravalley and intervalley scat¬ 
tering contributions to the relaxation time of the charge carriers, as these two valleys are 
occupied by the carriers for n-type SL- M 0 S 2 . Here the first order deformation potential is 
also shown in addition to the other scattering mechanisms. The longitudinal acoustic (LA) 
and transverse acoustic (TA) modes contribute to quasiclastic intravalley scattering process 
through the acoustic deformation potential (ADP). The contribution from TA is non-zero 
because of inclusion of Umklapp process and results in a coupling of highly anisotropic na¬ 
ture [69]. From Fig. ]24](a), it is clear that at low carrier energy the total scattering rate is 
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FIG. 23. The exponent ( 7 ) for mobility as a function of phonon energies (Hoj). Exponent is given 
by Hbuik = Aio(T/To) -7 . The graphs a, b and c correspond to the scattering by homopolar phonon 
modes, by polar phonon modes and by acoustic phonon modes, respectively. Here, bulk-MoS 2 
corresponds to huj=60 meV and 7 = 2.6 (indicated by an arrow along the vertical axis). Taken from 
Ref. [67]. 

dominated by ADP. At higher energy the total scattering rate is dominated by zero-order 
optical deformation potential (ODP) and the scattering via Frohlich interaction; the sudden 
jump in the optical scattering rates is because of onset of the emission of optical phonons 
by the charge carriers. At all energy scale, the first-order ODP scattering is one order of 
magnitude smaller than the other scattering rates. 

Similar to the bulk case, the mobility due to ADP is /x“£ oustlc oc r p-i acoustlc |ggj. w here 
7 acoustic _ -p 7 he calculated temperature dependent mobility for IL-M 0 S 2 is shown in 
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FIG. 24. (color online)-(a) The first-principles calculations of different electron-phonon scattering 
rates (and the total) at room temperature versus the charge carrier energy. Here ADP=Aa, Zero- 
order ODP=D^ and First-order ODP=D^. The dashed line shows the fitted deformation potential 
parameters given by = A \q for acoustic modes and M (l \ = D^q for all the other modes by 
choosing corresponding deformation potentials. The kinks in the curves for the optical phonon 
processes are attributed to the onset of phonon emissions, (b) Calculated mobility for different 
temperatures. The shaded area denotes the variations of the mobility due to the uncertainty (10%) 
in the calculated deformation potentials. For comparison, the mobility because of the presence of 
the acoustic modes only is plotted, where /x ~ 1/T. The red solid line represents the improvement 
of the mobility due to the quenching of the homopolar modes. Taken from Ref. [69 ]. 


Fig. 24 [b). At lower temperature (~ 100 K) the value of 7 = 1 . At higher temperatures the 
mobility follows according to fi°^ lcal oc T~ 7 , where 7 = 1.69 due to zero-order ODP and 
Frohlich scattering [69]. The calculated mobility at room temperature is of ~ 400 cm 2 /V 
sec. If we exclude the zero-order ODP (homopolar mode), the resulting mobility is increased 


by rs_x 70 cm 2 /V sec and the reduced 7=1.52, as shown in Fig. ]24[b) for the quenched case; 
although the effect of quenching on the mobility is not much. 


Comparison with the experimentally observed mobility: The experimentally 
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FIG. 25. (color online)-The temperature dependent Field effect mobility for single layer M 0 S 2 
deposited on Si 02 /Si substrate. The back gate voltage (Vbg) is in the range of 30-40 V to extract 
the mobility.. In the high-T range, the fitted value of the exponent ( 7 ) is 1.4. Taken from Ref. 

Ell- 


observed (at room temperature) low mobility in the range of 1-50 cm 2 /V sec [10] 170, El] 
reported earlier is due to the charge impurity scattering and increased upto 200 cm 2 /V sec 
[TO] by depositing a high dielectric material (HfCb) on top of the SL-M 0 S 2 . This enhance¬ 
ment of the mobility is partially attributed to the strong damping of Coulombic scattering 
from charge impurities [72]. The enhancement of mobility is also achieved by polymer 
(PEO+LiClCb) top gating from ~ 0.1 to ~ 150 cm 2 /V sec [73]. The difference between the 
calculated value (~ 410 cm 2 /V sec) and the experimentally observed (~ 200 cm 2 /V sec) 
is may be due to the defects and surface optical phonon scattering [69]. The temperature 


dependent mobility with the single-gate configuration is shown in Fig. 25. The field effect 
mobility is given by fj, — [da / dVb g ][L/WCb g }, where a denotes conductance, V ) >g applied 
back gate voltage, L channel length, W the width of the channel and Q )ff represents back 
gate capacitance associated with the Si02/Si layer. At ~ 200K, the mobility reaches a peak 
value of ~ 18 cm 2 /V sec. The decrease of mobility with temperature from 200K to 4K is 
attributed to the well known scattering from charged impurities. Above 200K, the mobility 
decreases with the increase in temperature. This high-T limited mobility due to phonons is 


34 






fitted according to jisL ~ T 7 , where 7 ~1.4. This is in good agreement with the predicted 
value of 7 «1.65 [69j. With the dual-gate configuration, the observed mobility is plotted 


against temperature, as shown in Fig. 26 The mobility varies from ~ 168 cm 2 /V sec at 
4K to ~ 60 cm 2 /V sec at 240K. The mobility with dual-gate configuration is quite different 
than that with the single-gate configuration, where the mobility monotonously decreases 
with decreasing temperature. This distinct feature with dual-gate device is attributed to 
the strong damping of charged impurity scattering in presence of dielectric media and to 
the metallic top gate which also changes the dielectric environment for the SL-M 0 S 2 [31] • In 
the high temperature (100-300K), the fitted value of mobility due to the phonon scattering 
shows 7 ~ 0.73. All the other dual-gate devices of SL-M 0 S 2 show the variations of 7 from 
0.3 to 0.75 [31]. These observed values are much less than the predicted 7 ~ 1.52 [69] . 
In addition to the quenching of homopolar mode due to the dielectric environment, other 
phonon scattering mechanisms are also reduced which needs further theoretical work to 
clarify this discrepancy. 
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FIG. 26. (color online)-The temperature dependent Field effect mobility (/i=[d< 7 /dV;, 9 ][L/WCig) 
for single layer M 0 S 2 deposited on Si 02 /Si substrate. Applied top gate voltage (Vtg) is in the 
range of 3-5 V to extract the mobility. In the high-T range, the fitted value of the exponent ( 7 ) is 
0.73. Taken from Ref. [31]. 
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6.3 . Substrate effect towards improving the mobility of M 0 S 2 : suppression of 





FIG. 27. (color online)-Schematic diagrams of the four-probe devices (contacts made of Ti/Al 
electrodes) (a) without PMMA and (b) with PMMA (300 nm) deposited on 300 nm Si02/Si 
substrate. The measured conductivity (<r) of four MoS 2 -devices having PMMA with different 
thicknesses (1.5 nm, 6.5 nm, 47 nm and 80 nm) (c) with linear scale and (d) with the semi-log 
scale. Taken from Ref. cu¬ 
lt has been established that the presence of the top gate material on M 0 S 2 enhances 
the mobility of the charge carriers [TOl [3T1 73j. Before putting the top-gate material the 
SL-M 0 S 2 deposited on bare Si0 2 /Si substrate shows a less amount of mobility in the range 
of 1-50 cm 2 /V sec [TUI 1701 71] . In order to investigate the substrate effect on the measured 
mobility, two different kind of devices are considered for the various thickness (1-80 nm) of 
M 0 S 2 : (i) the M 0 S 2 flakes are deposited on Si02/Si substrate (shown in Fig. [27fi) and, (ii) 
those deposited on polymethyl methacrylate (PMMA) (shown in Fig. |27jp). Fig. [27|(c) shows 
the observed room temperature conductivity (c) with the applied back gate voltage (V;, 9 ) 
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for four different thickness’s of sample deposited on 300 nm PMMA m■ Fig. pWd) shows 
the semi-logarithmic behavior of the conductivity. While the 6.5 and 47 nm thick M 0 S 2 
show the ambipolar behavior, 1.5 nm thick ( 2 L-M 0 S 2 ) shows unipolar nature with mobility 
30 cm 2 /V sec. The mobilities for the 6.5 nm thick sample are given by /z e (6.5) ~ 68 cm 2 /V 
sec and /x/ l (6.5) ~ 1 cm 2 /V sec for electrons and holes, respectively; for 47 nm thick M 0 S 2 , 
these are given by / u e (47) ~ 270 cm 2 /V sec and fih{ 47) ~ 480 cm 2 /V sec. It is clear from 


Fig. 27(d) that the region (off-state), separating electron and hole conductance, decreases 
with the increasing thickness of the sample; the off-state conductance also increases with 
the thickness. For 80 nm thick sample, there is no off-state indicating the dominance of 
conductance from the bulk m- 
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FIG. 28. (color online)-The field effect mobility (at 300K) versus thickness of the M 0 S 2 flakes for 
25 devices with PMMA and 6 devices without PMMA. For PMMA supported devices (range of \j, g 
± 150V), electron (solid squares) and hole (empty squares) mobilities are connected with dashed 
lines. For the Si02 supported devices (range of Vf )fJ ± 75V), only electron mobility (red circle) is 
shown. Taken from Ref. m- 


The thickness dependent mobility at room temperature for two different kind of device 


configurations is shown in Fig. [28} The ambipolar nature has been marked by a dashed line 
connecting the hole and electron mobility. From Fig. 


all the devices can be grouped in 


two categories: (a) unipolar type-shown by all the devices deposited on Si0 2 /Si substrate 
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and by all thin samples (1-30 nm) supported by PMMA ; (b) ambipolar type-observed for 
all the thick sample of width in the range of 40 nm to 70 nm deposited on PMMA dj. 
The observed mobility peaks are for ~ 50 nm thick M 0 S 2 samples of ambipolar type and 


are given by /W e (50) ~ 470 and ^(47) ~ 480 in cm 2 /V sec. From Fig. 28 it can be inferred 
that for most of the devices of ambipolar type, the [ih is larger than the /i e . Hence, the 
multilayer M 0 S 2 for the p-type operation can be chosen for faster electronics. The thick¬ 
ness independent mobility shown by unipolar-type devices deposited on Si0 2 /Si substrate 
is in the range of 30-60 cm 2 /V sec. The observed higher mobility at room temperature ( 
/i e (50) rv_/ 470 and /Th(47) ~ 480 in cm 2 /V sec) for the multilayer M 0 S 2 is larger than the 
theoretically predicted phonon limited mobility for SL-M 0 S 2 [69j and for the thick layer of 
MoS 2 (~ 100 cm 2 /V sec) [[75]. 
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FIG. 29. (color online)-The conductivity with back gate voltages (a) for PMMA supported device 
of 53 nm thick M 0 S 2 and (b) for SiCH supported device of 42 nm thick MoS 2 - The red and blue 
curves show the conductivity before and after PMMA coating (from top) on M 0 S 2 flakes. The 
insets show the schematics of the devices. Taken from Ref. ITU- 


In case of M 0 S 2 devices, the charge carriers are confined close to the dielectric interface 
(within few nm) used for the back gating [75] . The thickness dependent enhancement of the 
carriers mobility can be attributed to the role of the additional M 0 S 2 layers to screen the 
long range disorder m To check the screening effect, PMMA coating has been done onto 


the top of the devices for the two different types of configurations (see Figs. 27 1 and b). The 
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resulting measured conductivity is shown in Figs. [29|a) and (b). While the PMMA coating 
helps to screen the long range scattering for the PMMA supported devices, but the effect is 
little for SiOo/Si supported devices. The maximum improvement mobility observed with the 
PMMA coating is ~ more than 300 % for 60 nm thick sample 0- The observed variation 
of the mobility on the thickness or to the PMMA coating is attributed to the presence of 
short range disorder (rough surface) for the SiC^/Si supported devices [73]. 

6.4 . Reduction of Schottky barrier at metal-semiconductor contact to improve 
the mobility of M 0 S 2 



FIG. 30. (color online)-The drain-source current versus back gate voltage for 6 nm thick M 0 S 2 . 
The threshold voltage, Vt/i= -6.0, -1.0, 1.5 and 4.0 V for Sc, Ti, Ni and Pt contact devices, 
respectively. The inset shows the alignment of the Fermi levels of four metals with the conduction 
band edge (irrespective of different work functions of the metals) according to the observed n-type 
device characteristics. Taken from Ref. Ell¬ 
in all the previous studies, attention has been paid to the interface physics i.e. the 
interface between gate dielectric and the sample, the interface between the substrate and 
the sample to improve the mobility and hence device performance. However, the presence 
of Schottky barrier between metal contacts and semiconductor can affect the mobility sig¬ 
nificantly. In order to understand the Schottky barrier contacts and their role in improving 
the mobility, four different metals scandium (work function $^=3.5 eV), titanium ($^=4.3 
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eV), nickel (<3 >m= 5.0 eV) and platinum (<E>m= 5.9 eV) are used to make contacts with the 


shows the experimentally measured transfer characteristics for the 6 nm thick M 0 S 2 device 
m All the four devices show the unipolar n-type operation, irrespective of the different 
work functions of the four metals. On the positive side of Vgs-Vth, there is a clear decrease 
in on-current from Sc to Pt contacts, and this is consistent with the formation of Schottky 
barrier at the metal to channel interface. 


multilayer MoS 2 (electron affinity = 4.0 eV [7B] ) for the devices using a back gate. Fig. 31 


' r FB 'iH 



FIG. 31. (color online)-(a) The drain-source current versus back gate voltage at different tempera¬ 
tures showing three distinct regions of operations: below flat-band voltage (Vfb), below threshold 
voltage and above threshold voltage. The insets show the energy-band diagrams in the corre¬ 
sponding regions, (b) The drain-source current versus the inverse-temperature (1/T) (taken from 
a), (c) The effective barrier height (&b from b) as a function of gate voltage and the extraction of 
true Schottky barrier (&sb) for Ni-contact device, (d) All &sb versus work functions for all the 
different metal contact (Sc, Ti, Ni and Pt) devices. Taken from Ref. [TT j. 
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In order to estimate the true height of the Schottky barrier (&sb) quantitatively, the 
temperature dependence (100-400K) transfer characteristics is measured for the Ni contact 
as shown in Fig. |3lj(a). In the subthreshold region, the number of charge carriers are less and 
hence diffusive temperature dependent thermionic current ( Ithermionic ) and tunnelling current 
(I tunneling) become important. The temperature dependent current (Ids) is then plotted for 
various back gate voltages, as shown in Fig. |3lj(b). The thermionic emission current is fitted 
to extract the effective barrier height (<h B ) with the well know formula Ids— I thermionic 
= AT 2 exp[-^p}[ 1 — exp(jf ^)]; where, A is the Richardson’s constant, K B the Boltzmann 
constant, e the electronic charge. Fig.[3l|(c) shows the variation of the extracted & B with the 
applied Vgs- The effective barrier height = &sb — C e ff(Vcs — Vfb ) [Uj, for Vgs < Vfb 
and $d = <3>sb — A, for Vgs > Vfb- The factor C e ff = 1/(1 + Cit + Cch ) represents the 
band movement of the channel with Vgs, where C lt and C ox denotes the capacitance for 
the interface trap charges and oxide layers, respectively. A is a positive quantity but not 
constant. V fb is the flat band voltage. The Ithermionic is the only current through the channel 


for the negative Vgs < Vfb (see Fig. 31 1 ); and hence, the linear relationship between 
and as shown in Fig. [3l](c) . When Vgs > Vfb, the Itunneiing starts to flow through the 
channel and this current is not included in the fitted equation for Ithermionic- And because 
of this fact, the linear relationship between $d and &sb does not hold and starts to deviate 
after Vfb, as shown in Fig. [3l](c) . The extracted value of & S b is 0.15 eV fo Ni contact. It 
must be noteworthy that this is the tunneling current which is responsible for getting ohmic 
contact even for the higher work function metals including gold ($^=5.4 eV) [IT] . Fig.Qd) 
represents the measured &sb for different metals, which are 30, 50, 150 and 230 meV for Sc, 
Ti, Ni and Pt, respectively. The fitted dashed line gives the slope d^ss/d^M ~ 0.1 (0.27 for 
Si [77]) implying that there is a strong pining of Fermi level at the semiconductor interface 
before the metal contact is made. This pining of Fermi level near neutral level is related 
to the charge impurities present on the semiconductor-metal interface, which is usually the 
case with the covalent semiconductors m This explains clearly why only n-type device 
characteristics are observed irrespective of high work functions of the metals used for making 
contacts. 


To understand the effect of measured Schottky barrier on mobility, the output and 
transfer characteristics for 10 nm thin M 0 S 2 are shown in Fig. [32|(a). The exponential 
rise (see the inset) of I DS with V B s f° r the Pt contact is in accordance with the Schottky 
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(a) 



(b) 



FIG. 32. (color online)-(a) The drain-source current versus back gate voltage (T=300K) for Sc, Ti, 
Ni and Pt-contact devices. Here Vgs-Vth =4 V corresponds to ~ 1.4 x 10 12 /cm 2 . The value of 
the back gate capacitance is ~ 3 x 10 -4 F/m 2 . (b) The extracted held effect mobility for different 
thicknesses of the sample with Sc-contact devices. The dashed line shows the fitting of the data 
according to the resistor network model, where two model parameters are Thomas-Fermi screening 
length (Xth ~ 3 nm) and the interlayer conductivity {(Tint = 1 mS). The solid line represents the 
calculated mobility without taking The inset shows the schematic of the resistor network 

model. Taken from Ref. [ 113 - 

barrier of 230 meV. For Ni-contacted devices having $5^=150 meV, it shows the linearity. 
Since at 300K, the temperature assisted tunnelling current is responsible for this observed 
linear behavior between the current and drain voltage; the issue of getting ohmic con¬ 
tact is irrelevant m- The measured field-effect mobilities from transfer characteristics 
are 21, 36, 125 and 184 in cm 2 /V sec for the devices having Pt, Ni, Ti and Sc contacts, 
respectively HU- The measured low mobility M EH EH E3] (< 50 cm 2 /V sec) using Au or 
Ti/Au contacts without top gate material is in agreement with the effect of Schottky barrier. 


Fig. 32(b) shows the variations of the mobility with the thickness of the sample in the 
range of 2-70 nm for the device having Sc contacts. The maximum mobility observed 
corresponds to the ~ 10 nm thin M 0 S 2 device. The observed mobility is fitted with the 


resistor network model (see the inset of Fig. 32 d). For thicker samples, the access to the 
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FIG. 33. (color online)-The drain-source current versus drain-source voltage (for different top gate 
voltages) at T=300K for 10 nm thickness of the M 0 S 2 layer (Sc-contact). The improvement of the 
transport coefficients are achieved by depositing a top-gate material (AI 2 O 3 ) of 15 nm thickness. 
The inset shows the corresponding transfer characteristics for different drain-source voltages. Taken 
from Ref. [Hi- 


lower layer involves the interlayer resistance. For thinner samples, the screening of the 
substrate is much less than the thicker one. The observed peak is, therefore, related to the 
competition between the interlayer conductivity (<jj nt ) and the screening of the substrate 


mi- Fig. shows the the output and transfer characteristics for 10 nm thin M 0 S 2 with 
the addition of 15 nm dielectric material (AI 2 O 3 ) on the top of the sample. The held effect 
mobility due to the screening of long range scattering is now enhanced to the value of ~ 700 
cm 2 /V sec from 184 cm 2 /V sec for Sc contacted device. This huge increment is attributed 
to the (i) use of low Schottky barrier ( < h§£=30 meV), (ii) choosing of optimum value of layer 
thickness (10 nm) which corresponds to the optimization of suppression of substrate effect 
as well as the interlayer resistor, and (iii) screening of the long range Coulomb interaction. 
Finally, a high saturation current density ~ 240 [iAjfi m at the carrier density of ~ 2.8 x 
10 12 /cm 2 and the higher transconductance ~ 4.7 )iS//mi at Vds= 1-0 V are measured for a 5 
/mi channel length. In order to get a higher mobility (> 700 cm 2 /V sec), it is desirable to do 
experiment on a device having Sc contact deposited on a high dielectric material instead of 
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Si02/Si substrate (may be PMMA |74j ) with high dielectric top gating material by varying 
only the thickness of the sample. 


7. PHONON-RENORMALIZATION IN A SINGLE LAYER MOS 2 TRANSISTOR: 
ASPECTS OF SYMMETRY 

It has been established that the high temperature (~ 200 to 300 K) phonon limited mo¬ 
bility in single and bulk-MoS 2 (homopolar mode A\ g ~ 408 cm -1 ) is because of dominant 
scattering of charge carriers by the optical phonons HH E7j. Electron-phonon scattering 
not only affects carrier momentum relaxation, but also affects the frequency and the full 
width at half maximum (FWHM) of phonons via renormalization of their self energy; the 
real part of the self energy is related to the frequency shift and imaginary part to the 
FWHM [EH]. The non-destructive Raman spectroscopy is exploited extensively to probe the 
electron-phonon coupling (EPC) quantitatively for single layer graphene [79H8T] and bi-layer 
graphene [821 - 184] . To quantify the EPC for single layer MoS 2 , in-situ gate voltage dependent 
Raman experiments have been reported [7S]. The schematic diagram for the device is shown 
in Fig. |34|(a). A mixture of LiClCR and polyethylene oxide (PEO) in 1:8 weight ratio is 
used as a top gating material. Fig. [34](b) shows the atomic force microscopy (AFM) image 
in contact mode to measure the height (0.7 nm) of the MoS 2 . The output characteristics 
of the device is shown in Fig. |34](c). The observed slight non-linearity between I^s and 
V ns is due to the Schottky barrier for Au contacts m (as discussed earlier). The transfer 
characteristics in semi-log scale is shown in Fig. [34](d) . The n-type operation was attributed 
to the Fermi level pinning at the semiconductor interface due to interface charges HU • An 
extracted mobility was ~ 50 cm 2 /V sec and the on-off ratio was 10 5 (maximum observed ~ 

10 8 [ID]). 


The observed Raman spectra with different top gate voltages is shown in Fig. |35[a). In 
the spectral range of 350- 450 cm -1 , two Raman modes are observed. From Table-I, these 
two mode are A) ( rv_/ 402 cm x ) and (~ 382 cm x ) corresponding to the bulk A lfl and 
El g modes, respectively. We will continue to use the notation Ai 9 and E\ g for the observed 
two modes. The frequency variations of the two phonon modes with the top gating (Vtg) 

b). The ki g mode softens by 4 cm -1 and the E\ g mode shows almost 


are shown in Fig. 
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FIG. 34. (color online)-(a) Schematic of the single layer M 0 S 2 device (50 nm thick Au contact 
used). Here a;,; and uj s represent the incoming and scattered photons, respectively. The length of 
the device is ~ 2.5 /jm and width is ~ 1.5 /xm. (b) The AFM height profile of the single layer. 
The right inset (AFM image) showing the dotted line on the sample along which the AFM height 
profile has been taken, (c) The output characteristics of the device at different top gate voltages, 
(d) The transfer characteristics in semi-log scale (the inset showing in linear scale) at Vd5=200 
mV. It shows the on-off ratio ~ 10 5 . Taken from Ref. [75] . 


no change (~ 0.6 cm 1 ) in frequency upto the maximum doping of electron concentration 


1.8 x 10 13 /cm 2 (Vtg ~ 2 V). The FWHM of Ai g mode increases by 6 cm -1 (see Fig. 35:), 
but the FWHM of mode remains constant [78]. Similar to the detection principle of 
layer numbers of MoS 2 by measuring the uja,„ -w e i (46], the Accan be used to estimate 

9 2g 9 


the doping concentration of the MoS 2 flake. 


In order to understand the renormalization of the Ai g phonon mode, first-principles 
density functional theoretical (DFT) calculations have been carried out [78J. Single layer 
MoS 2 has a direct band gap of 1.9 eV at the K-point. The hybridized conduction band 
edge consists of dominated d 2 2 -type Mo states and the hybridized valence band edge states 
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FIG. 35. (color online)-(a) Raman spectra along with Lorentzian fits at different top gate voltages 
for the single layer MoS 2 - The changes in frequency and FWHM for the Ai s (b and d) and (c 
and e) modes as a function of top gate voltages. Taken from Ref. M- 

having d^-type [78], as shown in Figs. |36)(e) and (f). The square of the electronic wave- 

function (|'0( r )| 2 ) or the charge density near the valence band top and the conduction 

band bottom has the same full symmetry as the single layer M 0 S 2 crystal [78]. The 

EPC for a phonon mode v with momentum q and frequency is given by [85] \ qv = 

hjjqyN(cf) EkEqlCV x 5(ek+ q ,i — e/)x S(e kj — £f), such that the EPC matrix is given 

1 

by (fit+q.k = ( 2itfa q „ ) 2 (V’k + q,i | A V qu |^kj), where N(e/) is the density of states at the Fermi 
energy for electrons, e^j the electronic energy with momentum k in a band j, corre¬ 
sponding electronic wave function. Here AF qi/ denotes the change in the potential due to 
the vibrations of the lattice associated with a phonon (o; qi/ ). In the A lff phonon mode, the 
atoms vibrate in such a way that the total symmetry of the single layer M 0 S 2 does not 
change, and hence it corresponds to the identity representation. For n-type device opera¬ 
tions, the electrons are occupying gradually the conduction band edge states and for those 
electrons |'0(r)| 2 transforms according to the identity representation. Hence, EPC matrix 
element is non-zero, in which the product of AV] 4 lg and |-0(r)| J is involved. The irreducible 
representation for the E\ g mode is orthogonal to the A lff mode representation and hence the 
EPC matrix vanishes leading to no response of this mode to the channel doping ra- 

To understand dependence of FWHM on doping, the total FWHM is expressed as 
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FIG. 36. (color online)-Experimentally observed and the DFT calculated changes in frequency 
of the (a) Ai 9 and (b) E \ g modes versus electron concentrations (n). (c) Calculated EPC for the 
Ai g and E\ g modes, (d) Calculated two ratios for EPC [EPC(n / 0)/EPC(n=0)] and FWHM 
[FWHM totoi (n / 0)/FWHM totoi (n=0)] of the Ai fl mode. The charge densities (|V ; ( r )| 2 ) for the 
(e) top of the valence band and (b) bottom of the conduction band. Here, The dark gray (red) 
and black (blue) spheres indicate Mo and sulfur atoms, respectively. Taken from Ref. [75] . 


FWHM ioioi =FWHM E PC +FWHM an ; where FWHM an is due to anharmonic effects. The 
two ratios FWHM totai (n ^ 0)/FWHM totaZ (n=0) and EPC(n ^ 0)/EPC(n=0) have been 
calculated for the A lg mode and they show the same trend as shown in Figs. |36|(d). This 
explains that only EPC is responsible to broaden the FWHM of A l£ , mode with the doping. 
It is noteworthy that the top gating by an amount of ~ 1.8 x 10 13 /cm 2 leads to hardening of 
the G phonon mode by ~ 10 cm -1 in case of graphene ra. while the same condition gives rise 
to softening (by ~ 4 cnr 1 ) of Ai 9 mode for single layer M 0 S 2 . In case of graphene adiabatic 
approximation fails, whereas in case of M 0 S 2 , EPC is within the adiabatic approximation. 
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8. RESONANT RAMAN SCATTERING OF BULK MOS 2 TUNED BY TEMPER¬ 
ATURE AND PRESSURE: STABILITY AND CONTRIBUTIONS OF A AND B 
EXCITONS TO RESONANCE 


Resonant Raman scattering at ambient conditions on single and multilayer MoS 2 has 
been discussed m, where 633 nm red laser line was used for excitation. It has been shown 
that there exists two excitons corresponding to direct transitions at the K-point |64j . The 
exciton A having binding energy 42 meV corresponds to K 4 —» K, 5 optical transition, whereas 
B exciton of binding energy 134 meV corresponding to —> K 5 optical transition. The en¬ 
ergy difference between the A and B excitons is due to the spin-orbit splitting and interlayer 
interaction H- Since, the pressure and temperature variations affect the lifetime of the 
excitons, we will address their stability and the limits of applied pressure and temperature 
to which the resonance effects can be observed. The appearance of all the Raman active 
modes and inactive modes (E 4u and B lu ) due to resonance effects is already discussed [47] hr 
the section V. Here, our focus will be on the V mode attributed to the two-phonon Raman 


process (see Fig. 20) such that it corresponds to the emission of a dispersive quasi-acoustic 


(QA) longitudinal phonon followed by emission of a dispersionless transverse optical mode 
(here c-band) [60]. The more interesting feature about the b mode is that the observed 
frequency for the Stokes (cuf) and anti-Stokes (co^ s ) spectra is not the same and they differ 
by ~ 4.7 cm -1 , as shown in Fig. 37 



FIG. 37. (color online)-(a) Temperature and (b) pressure dependent oj^ and u£- s as a function of 
the energy difference A E = El — Eai ■ Solid line in (a) indicates the fitting for . Taken from 
Ref. [SB]. 
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Since, longitudinal phonons are strongly coupled to excitons than the transverse ones, the 
contribution of QA phonons will be more dominant to the b mode [60]. The rate at which 
frequency decreases with temperature ( du/dT ) for the b mode is ~ 2.66 x 10 -2 , whereas 
for all the other modes it is in the range of ~ 1.0-1.6 x 10 -2 . Fig. [37|(a) and (b) shows 
the comparison between the shifting of the peak position of the b band for the Stokes and 
anti-Stokes scattering as a function of temperature and pressure, respectively. Here, the uj^ s 
is redshifted compared to the cuf part. 

Since, the intermediate exciton state plays an important role for the observed b band 
and from the above experimental findings, three major things need to be addressed: (i) 
identification of tuned excitonic states with respect to temperature and pressure to observe 
resonance effects for Stokes as well as anti-Stokes processes, (ii) the stability of the A and 
B excitons and, (iii) the observed redshifted frequency ~ 4.7 cm -1 . As we already 

noted [78], the d 2 2 ( this orbital points along the c-axis) character of the K, 5 final state 
predominantly couples the A lff phonon with the excitons than the mode. This means 
that the Ai 9 phonon mode is much more sensitive to the resonance effect than the E\ g . Hence 
the resonance effect will be captured from the intensity ratio by comparing [I a U JI e) 2 ] reso 
with the [Ia^JIeI ] non ~ reso (i.e. with 514 nm laser line). Here, the measured [Ia u J I e 1 2 } reso 
implies the Raman scattering cross-section. The intensity ratio [J^/Igi } reso is denoted by 
R 5 for Stokes process and by R: 45 for anti-Stokes process. 

The probability of the Raman scattering of a phonon by an excitonic Is state is given by 
PH 186] 

(0| H eR {u L ) |1) (1| H e _ ion |1) (1| H eR {u a ) |0) 2 

Here, 77 e # and H e _ ion represents the electron-photon and electron-ion interaction Hamil¬ 
tonian, respectively. The states indicated by 0 and 1 correspond to the ground and the 
exciton states (Ai and Bi), respectively. The incident and scattered laser frequencies are 
given by col and co s , respectively. E,; = Hul or =hu s depending on incoming or outgoing 
resonance, respectively. It has been shown experimentally that the Raman scattering cross 
section for the b-band shows an peak at E^, + hujf, (for the A exciton), but does not cor¬ 
respond to energy of the exciton Is state i.e. at E^ [60]. Here, the outgoing probability 
for the Stokes and anti-Stokes Raman scattering of Ai 9 mode is given by Pf and 
corresponding to Ei + fvujA lg and Ei — hujA lg energy levels, respectively. Ho et al [88] have 


Pphonon ~ (2"7r/i) 
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measured the temperature dependence of the exciton energies by fitting a formula given by 
Ei(T ) = E iPo — dip[ 1 + j T ^—\ ]; where, d iP and Q iP represent the exciton-phonon in¬ 

teraction strength and the average phonon temperature, respectively. The measured values 
of E iPo , di P and 0jp are for Ai are 1.976 eV, 46 meV and 220 K, respectively; and for Bi 
those are 2.179 eV, 42 meV and 200 K, respectively [ 88 ] for T < 300K. The broadening 
parameter has also been measured by fitting the equation T(T) = T 0 + e^e^/T)-! > where 
Ylo — 75 meV, represents the exciton-LO phonon interaction strength and ©lo= 560 K. 
The T 0 for Ai and Bi excitons are 18.0 and 37.4 meV, respectively. Hence, it is clear that 
the broadening is because of exciton-optical phonon interactions [ 88 ] . 


The temperature dependent binding energies of the Ai and Bi excitons are shown in 
Fig. 38'a). The measured intensity ratio R s and R: 45 , and the calculated probabilities Pf 
and are also plotted (blue and red solid lines). The outgoing resonance energies for 
the Ai g mode are also shown for Stokes and anti-Stokes scattering processes of Ai and Bi 
excitons. The interactions of excitons with the acoustic and optical phonons are responsible 
for the observed temperature dependent shift; and also for the variations of the lattice 
parameters. Since the exciton-optical phonon interaction holds for the high temperature, the 
fitted parameters are extended to 550 K. The agreement between the theory and experiment 
is reasonable. From Fig. |38](a), it is clear that at T=0K the excited laser line (Ep=1.96 eV) 
lies between Eai and Eai + hujAig- As we increase the temperature the Eai + fr^Aig (Stokes- 
process) approaches to and around ~ 95 K, outgoing resonance shows an peak for the 
R 5 i.e. for the Raman cross section. Now with the further increase of temperature for 
T > 175 K (marked as a on the upper scale) Eai + kiOAig goes away from the E L and 
E B \ — hujAig approaches E L and continues upto ~ 260 K (marked as (3). In the boundary of 
175/1 <T< 260/1, the resonance effects are of mixed type: due to the proximity effects of 
El to the Stokes/anti-Stokes shifted outgoing resonance for the Ai and Bi excitons. With 
increase of T further, R: s goes to a minimum value around 350K, whereas R : 45 reaches a 
maximum value at 400K (marked as 7 ). For T > 520 K (marked as <5), R/ s and R ^ 5 
follows the non-resonant temperature dependence behavior of exp(—hujAig/kT) [ 86 ]. 

Connell et al has shown experimentally that the binding energy of the Ai exciton decreases 
with pressure at the rate of ~ 7 meV/GPa at room temperature [89]. The decreasing of 
the exciton binding energy is associated with the pressure induced free charge carriers and 
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(a) (b) 



FIG. 38. (color online)-(a) (corresponding to the left scale) Temperature (measured for T < 300K 
and then extended to 550K) dependent binding energies (E^, and Es x ) of the Ai and Bi excitons. 
(b) Pressure dependent binding energy for the Ai exciton is plotted (corresponding to the left 
scale). Here AS and S represent the outgoing resonance energies E) + hojA lg and Ej — huoA lg for 
anti-Stokes and Stokes Raman scattering, respectively. The measured intensity ratios R 5 for Stokes 
process, R" 45 for anti-Stokes process and the corresponding calculated probabilities ( P s and P AS ) 
are plotted together (corresponding to the right scale) in (a) and (b). The inset in (b) shows the 
R 5 for pressures above 10 GPa. The upper scale in (a) showing a-5 represent different temperature 
regions of operations discussed in the text. Taken from Ref. [86]. 

hence more screening. Therefore, at ~ 6 GPa and higher values, the Ai exciton cease to 
exist. The pressure variations of the binding energy (~ 0.02 eV/GPa [86]), Stokes and anti- 
Stokes shifted energies of the Ai exciton are plotted in Fig. |38|(b) along with the measured 
(R 5 and R A5 ) and calculated (Pf and P A ^) resonance Raman scattering cross sections. 
The pressure dependent broadening parameter T(P) ~ 35 + 3 P (in meV) has been taken 
to calculate the P f and P A ^ values. According to calculations, the peak positions for 
the P f and P A ^ g are ~ 3.8 GPa and ~ 6.3 GPa, respectively [HE]- Since, Eai + huiAig is 
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closer to at low pressure regime, R s increases and peaks at ~ 3.8 GPa in accordance 
with the theoretical value. At higher pressure, E B1 — hu>Ai g approaches to (not shown in 
the figure) and hence, R/' s increases to the peak value at ~ 6 GPa. The deviation of P A ^ 
from R as for the pressure > 6 GPa is because of the dissociation of the Ai exciton into 
free electron-hole pairs [BE]. Although the Ai exciton is unstable in the high pressure range 
(> 6 GPa), the measured R 5 values with the increased pressure (P> 10 GPa) is shown in 


the inset of Fig. 38 b). Livneh et al attributed this observed high pressure regime of R s 
to the pressure induced change of electronic band structure, and by that the enhancement 
of the resonant sensitive A ls mode couplings. This point needs to be cleared further from 
theoretical point of view. 

To address the stability of the B exciton, Livneh et al [86] compared the observed intensity 
ratio I A 7h s ' and the calculated ratio Pai 9 / Pai 9 (normalized by n(uj p h) /[n(uj p h) + 1], where 
n(cOph) = l/[exp(cu p h ) — 1] is the Bose-Einstcin factor) of the Ai 9 and E\ g modes for anti- 
Stokes and Stokes Raman scattering as a function of 1/T and proposed that as long as the 
ratio 1 AS /I s of the Ai s mode follows the resonant feature, Bi exciton is stable; and hence it 
is unstable for T > 520 K. 

It has been observed that the b-band shifts to lower frequency with the increase of the laser 
energy and reaches to the frequency of A\ g mode [60J. This dispersive mode was assigned as 
the two phonon process having dominant contribution from the dispersive QA mode. In order 
to address the observed co^-lo as ~ 4.7 cm -1 feature, the temperature and pressure dependent 
cuf and u AS are plotted with the energy difference A E = E L — Eai in Figs. |37[a) and 


(b), respectively. The temperature dependent through A E data is fitted (see Fig. 37 1 ) 
according to the equation hub(q z , T) = A E{T) — h 2 (|g z | ± c<^Lyf ) 2 l^^-u + ^to{T )• Here, 
no represents the ordinary refractive index. M fj represents the effective mass of the A\ 


exciton along the c-axis. The ± sign refers to the a and f3 processes (see Fig. 20). The 
q z and k, represent wavevector of the QA mode and the exciton pointing along the c-axis, 
respectively. The value of fau 6 (0, 300) = 52.4 cm -1 is used as a single fitted parameter. The 
extracted value of the frequency for the QA mode ~ 56 cm -1 is consistent with the earlier 
observations from neutron scattering [[63]. The calculated temperature coefficient of the QA 


mode is given by y « -2 x 10 3 in cm 3 /K for A E < 110 meV [86]; where, 

Uqa = Mb — uj c . The observed ) P _ 0 ~ —2.66 x 10 -2 in cm _1 /K is much larger than 


the 


dui, 


QA 


dT 


P =0 


Hence, the observed high temperature coefficient of the b-band compared 
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to all the observed mode is attributed to the tuning of the exciton energy with respect 
to the temperature [HE]. The largest shift between the Stokes and anti-Stokes scattering 


of the b-band for which uj^=uj^ s (see Fig. 37 1 ) is given by 50 meV for A E < 110 meV, 

J. Finally, is 


which is close to the energy of c-band (TO mode ~ 377 cm -1 ] 
attributed to the reverse order processes for the Stokes and ant-Stokes scattering: the inner 
photon-like state scatters to the outer exciton-polariton state by absorbing a TO phonon 
and then scatters back to inner high-energy photon-like state by absorbing a QA mode, i.e. 
for anti-Stokes process it is TO+QA order, whereas for Stokes scattering it is of QA+TO 
order. So, the redshift between cuf and co^ s is due to the order of the involvement of the 
dominant quasi-acoustic (QA) phonon in resonant two-phonon scattering process. 


9. PHOTOLUMINESCENCE OF MOS 2 : CONTROL OF VALLEY POLARIZA¬ 
TION BY OPTICAL HELICITY AND OBSERVATION OF TRIONS 

9.1 . Layer Dependence 


The phenomena of radiative recombination of optically excited electron (e) and hole (h) 
pairs of a band-gap semiconductor is known as photolumincscence (PL). PL study has been 
employed extensively to probe excitons and the band gap of a semiconductor. As discussed 
that the bulk MoS 2 undergoes a transition from an indirect band gap to a direct one in the 
monolayer limit [27], PL experiments are useful to probe the transition with the number of 
layers. Fig. |39|(a) shows the comparison of the PL spectra for the monolayer and bilayer 
MoS 2 . The striking difference with the addition of another layer is observed for the PL 
spectra, which is characterized by the photoluminescence quantum yield (QY). The PL 


spectra with the layer N=1 to 6 is shown in Fig. 39 [b). The observed peaks are marked as 


’A’, ’B’ and ’I’. The extracted QY for the bright PL of the monolayer is ~ 4 x 1CT 3 ; whereas 
for N=2-6, it is of ~ 10 ” 5 to 1CT 6 [27]. The inset of Fig. |39](a) shows the variations of the 
QY with the layer numbers. The observed peak A centered at ~ 1.90 eV corresponds to the 
direct exciton transition at the K-point of the BZ (K 4 —>■ K 5 optical transition); while the 
B peak centered at ~ 2.05 eV corresponds to the direct exciton transition Ki —> K 5 . The 
B peak is also observed for the monolayer MoS 2 deposited on Si0 2 /Si substrate [17| [39]. 
The A peak having width ~ 50 meV (for N=l) redshifts very little and broadens with the 
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thickness of the sample. The peak marked as I starts to appear with the bilayer and shifts 
to lower energy with the N. The broad peak 1 is attributed to phonon assisted PL peak of 
the indirect band gap MoS 2 . Since, the states near conduction band edge at the K-point 
of the BZ have d z 2 character and Mo atoms reside in between two planes of S atoms, these 
states are less affected for the interlayer couplings; whereas at the T point, both valence 
and conduction bands are linear combination of d-states (Mo) and p-states (S) and hence 
are more affected with the layer stacking. Therefore, indirect band gap is increased due to 
the confinement effects, direct band gap at the K-point remains almost constant. The peak 
energy of the A peak for monolayer and the peak energy of the I peak for all the other layers 
are plotted together in Fig. |39](c), which shows that the band gap decreases and approaches 
to the bulk value with the layers. 




FIG. 39. (color online)-(a) Room temperature PL spectra for suspended monolayer and bilayer 
samples. The inset shows the changes in QY with different layers (N=l-6). (b) The normalized 
PL spectra (with respect to the A-peak intensity) for all the layers of MoS 2 . Here, x 3 and X 10 
indicate the magnification of the I-peak in the corresponding spectra, (c) The A-peak energy for 
N=1 and the I-peak energy for N=2-6 plotted together showing the band-gap variations with the 
number of layers. The dashed line represents the indirect band gap of the bulk-MoS 2 . Taken from 
Ref. I2Zj. 

The direct band gap at K-point does not change with the layers, but the direct ex- 
citonic transition at the K-point gives different QY for the monolayer and the thick 
layer/bulk. In order to address this difference, the quantum efficiency for the MoS 2 is 
given by 77 « T rad /(T rad + r defect + T rdax ) [39]; where, T rad represents the radiative recom- 
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bination rate, Tdefect the rate of defect induced scattering, and T re i ax the rate of intraband 
relaxation of charge carriers. The phonon-assisted r re i ax is very large for the indirect band 
gap semiconductor. Because of the presence of the direct band gap of the monolayer MoS 2 , 
the decay rate via phonons decreases enormously i.e. r re i ax ~ 0. And r ra d is nearly constant 
for both the bulk and monolayer. Hence, the enhanced QY limited by Tdefect. for monolayer 
is due to the quenched channel of T re i ax with respect to the bulk value. 



1.6 1.8 2.0 2.2 2.4 

Energy (eV) 




Energy (eV) 


FIG. 40. (color online)-(a) The normalized PL spectra measured at 4K on monolayer MoS 2 
deposited on Si0 2 /Si substrate for two different laser (532 nm) powers. The two peaks are fitted 
with Gaussian functions, (b) The total PL intensity (integrated area of the two peaks) as a function 
of laser power, (c) The normalized PL spectra for different temperatures, (d) The energy of the 
H-peak versus temperature. Solid line represents the fitting of the Varshni formula. Taken from 
Ref. [90]. 
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FIG. 41. (color online)-(a) Time-resolved PL for different temperatures using 402 nm laser line, 
(b) The extracted decay rate and (c) the amplitude of the LLC. Taken from Ref. [901. 

9.2 . Temperature Dependence 


At the lowest temperature (~ 4 K), the observed PL spectra of monolayer M 0 S 2 for the 
two different laser powers (200 gW and 40 mW) are shown in Fig. |40|(a). The lowest broad 
peak is marked as ’L’ and the higher narrow peak is marked as ’H’. The L peak is separated 
from H by ~ 90 meV (A E s ). With the increasing power of the laser line by more than two 


orders, the shape of the spectra does not change. Fig. 40 b) shows the variation of the total 


integrated intensity of the two peaks with the laser power, which shows linear relationship. 
It does not show the onset of the PL emission due to the presence of the non-radiative 
recombination channel of charge carriers to the defects at the low power. The constant 
relative intensity of the two peaks with the laser power indicates that the two radiative 
recombination channel do not interfere each other with the increasing population of the 
charge carriers [90]. PL spectra at various temperatures (T) are shown in Fig. |40|(c). The 
two peaks redshift and broaden with the increased T. The intensity of the L peak decreases 
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and almost vanishes above 120K. The A E s remains constant before the vanishing of L peak. 
The redshift of H peak is plotted with T as shown in Fig. |40|d). The observed data is fitted 
with the Varshni formula, E g (T) = E g ( 0) — [aT 2 /(T + /?)]; where E ff (0)=1.874 eV, a— 5.9 x 
10 ~ 4 eV/K and {3= 430K. The H peak is attributed (tentatively) to the direct free excitonic 
direct transition at the K-point. The L peak is assigned due to the bound excitons with the 
defects. Due to the attraction of the neutral impurities/defects the binding energy of the 
bound-exciton decreases [ 68 ] . Because of the distribution of various defects with different 
energies in the sample, the observed L peak is broadened. The non-radiative recombination 
process increases with T and it is the responsible for the disassociation of the bound-exciton 
above 120K rather than the thermal activation processes (the binding energy is larger than 
the 26 meV). 

Fig. [4l|(a) shows the time resolved PL spectra at different temperatures. It is clear that 
the PL decays in 5 ps from 4K to 150K. This fast decay is in accordance with the observa¬ 
tion of relative intensity dependence on laser power at low T i.e. there is no non-radiative 
relaxation process of the optically generated e-h pairs down to the minimum energy. We 
can see from Fig. |4l|(a) that above 150K, there is a an appearance of another long-lived 
component (LLC). The extracted decay rate and the amplitude of the LLC are plotted in 


Figs. 41 b) and (c), respectively. The decay rate of LLC increase from 50 ps at 150K to 100 
ps at 270K and then decreases to 70 ps st 300K. The amplitude of LLC also increases with 
T, shown in Figs. |4l)(c). The appearance of LLC and its increased decay rate are attributed 
to the onset of non-radiative phonon assisted scattering channel of exciton-polaritons [90] . 


Plechinger et al [91J made two kind of samples of monolayer of M 0 S 2 to record the PL 
spectra: first kind (A and E) is of by mechanical exfoliation method on SiC^/Si substrate 
and the second type of sample (B and F) is covered with the oxide layer (Hf 02 or AI 2 O 3 ) 
of 15 nm thickness. The observed PL spectra at 4K, shown in Fig. |42|(a), present two kind 
of features: (1) A and E samples not covered with the oxide layer show two prominent 
peaks, and (2) B and F covered with the oxide layers show only one prominent peak. Since, 
the top oxide layers and the high temperatures used to make them remove most of the 
impurities from the surface, the low energy broad peak is attributed to the bound-excitons 
(to the surface impurities). The vertical dashed line in Fig. |42)(a) indicates the redshift of 
F sample covered with Hf 02 (e r ~ 20) and blueshift of B sample covered with AI 2 O 3 (e r 


57 





Energy (eV) Temp (K) 

FIG. 42. (color online)-PL spectra of (a) four samples (measured at 4K) and (b) three sample 
(measured at 240K). The A and E samples (in a) are not covered with oxide layers, whereas B 
and F are covered, (c) The high-energy peaks for three samples (shown in b) as a function of 
temperature. Taken from Ref. m- 


~ 10) compared to the without oxide layer samples (A and E). The screening of the long 
range Coulomb interaction does not play the role for the shifting of the B and F samples, 
otherwise for F sample the blueshift conld have been be more than the B sample. Fig. |42|(c) 
shows the temperature variations (4 to 240K) of the high energy peak for the two kind of 
samples. The largest shift (« 55 meV) occurs for the A and E type samples, whereas for 
the B sample it is of 34 meV and for F type it is of ~ 36 meV. Around 240K, the oxide 
covered samples show the blueshited peak compared to the without covered samples, shown 


in Fig. 42 b). Since, the atomic layer deposition process of making oxide layers involves 
temperature variations from 520K to below room temperature and the thermal expansion 
coefficient of M 0 S 2 is larger than the oxide layers, the sandwiched monolayer M 0 S 2 remain 
strained and hence the larger peak shifts are observed for the as-prepared samples than the 
oxide coated 
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9.3 . Effect of circular polarization of the incident light on photoluminescence: 


valley selectivity 


Even layers of M 0 S 2 has inversion symmetry, whereas odd layers and hence monolayer 
M 0 S 2 (D 37 J has no inversion symmetry. Because of no inversion symmetry, spin-orbit cou¬ 
pling splits the valence band of monolayer M 0 S 2 by ~ 160 meV at the K and K* points 
[92[ 193] (see Fig. |43)j). At the band edge of K points, the symmetry adapted hybridized 
d-states on Mo are given by \(p c ) = \d z 2 ) (1=0) and 1 <f) x ) = ~^(\d x 2 _ y 2 ) +ir\d xy )) (1=± 2). 
The hybridized p-states on S atom are \4> x ) = ~^{\Px) + \Py )) (1=± 1) ; where c and v 
represent the conduction and valence band, respectively. Here, r = ±1 indicates the valley 
index for K and K) points connected by time reversal symmetry (TRS). Those hybridized 
states interact with each other to form symmetry adapted linearly combined (SALC) states 
for the valence and conduction bands of monolayer M 0 S 2 . Now, the effect of the three¬ 
fold rotational symmetry (C 3 ) on those states are given by C 3 \v(K, K')) = \v{K,K')) 
and C 3 \c(K, K')) = e ^ 27 "/ 3 |c(A', K')) [32]• Therefore, the SALC states near the top of 
the conduction band dominated by d z 2 states on Mo carry an overall azimuthal quantum 
number m± = ±1 at K and K, whereas for the valence band SALC states it is given by 
rn± = 0. The optical selection rule for the direct transition at the K and K* becomes 
A m± = ±1. Since, the carriers have well defined angular momentum associated with the 
K and K' states, we need to have photons with particular hclicity (i.e circularly polarized 
light) to excite them. This phenomena of absorbing left-handed and right-handed circularly 
polarized lights by the two valleys is known as circular dichrosim (CD), see Fig. |43|(c). The 
TRS requires that Fq(k) = (-k) . Hence, we have the inherently coupled valley and the 

spins of the carriers because of the spin-filtered valence band splitting together with the 


TRS, as shown in Fig. 43(b). Fig. 43 (g) shows the PL spectrum (not polarization resolved) 
observed at T=14K using 2.33 eV laser excitation. The prominent feature around 1.9 eV is 
due to the A exciton complexes (AEC) consisting of two peaks: at higher energy it is due 
to the neutral exciton, and the lower component redshifted by 40 meV is attributed to the 
charged exciton (discussed later in details). The PL peak ~ 2.1 eV is due to the B exciton 
and ~ 1.8 eV arises because of excitons bound to the defects, as discussed earlier. 


I P cv (k) 1 2 — I P cv (k) 1 2 

The degree of circular polarization for the PL emission is given by r/(k, cu cv ) = |pA( k) 2 + | P A( k) _ ■ 
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FIG. 43. (color online)-(a) The hexagonal lattice structure for single layer M 0 S 2 and (d) the Bernal 
stacking for bi-layer M 0 S 2 • The valence and the conduction bands (marked by the corresponding 
total angular momentum) across the energy gap at the K and K points for (b) single layer and 
(e) for bilayer. Optical selection rules with the applied left-handed and right-handed circularly 
polarized lights (c) for single layer M 0 S 2 at the K and 1\ points, respectively; whereas for left- 
handed circularly polarized lights (f) in case of bilayer, (g) The unpolarized PL spectrum observed 
at 14K using 2.33 eV laser line for single layer MoS 2 - Taken from Ref. f29j. (h) The top and the 
bottom figures represent the conduction (pink) and the valence (blue) bands; and the middle one 
represents the calculated circular polarization r/(k, ui cv ) at different points of the Brillouin zone. 
The right scale bar represents the corresponding colors of the calculated ?y(k, uj cv ). Taken from 

Ref. | 32 ]. 

P±( k) representing the coupling strength with the circularly polarized lights (cr±) is given 
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by, P±( k) = k) ± Py V (k)}, where P™(k) = (\l/ c k| P" |\h„k). The calculated 


r/(k,Lu cv ) is shown in Fig. 43 h) indicating the exact value of ±1 at the K and K' points 


J. The polarization resolved (<r + and cr_) PL spectrum observed at T=14K using the 
left-handed circularly polarized (er_) lights (1.96 eV) is shown in Fig. |44)(a). The observed 
spectrum consists of only two peaks: The AEC peak is cr_ polarized and the peak due 
to bound-exciton is almost unpolarized. Mak et al ESI used a parameter to quantify the 
experimentally observed PL emission polarization, which is given by r] exp = 
where I(cr±) denotes the polarization resolved measured PL intensity. The observed value 
is given by 77 ea:p =1.00±0.05 and drops to ~ 0.05 below the photon energy 1.8 eV, as shown 


in Fig. 44 b). With the excitation of right-handed circularly polarized light (<7+), Mak et 


al observed the 7/ exp =-1.00±0.05 (not shown). In general, the quantity rj exp depends on the 
two timescales: (i) exciton lifetime given by r _1 = r r _1 + where r and nr represent the 
radiative and non-radiative recombination processes, respectively, and (ii) valley lifetime 
or hole-spin lifetime, t v [32]. hr terms of these two timescales, the helicity for the AEC is 


given by rf^ v = 1+ .,_ i ^ r t for A-neutral exciton and r/“ 


exp _ 


yj ex P 


1+2t a -/ t a-, 


for A-charged exciton 


J. Mak et al estimated that the observed helicity 77 e:Ep =1.00±0.05 is in accordance with 
the exciton lifetime r > 50 ps from the observed high QY on BN substrate and the valley 
lifetime r v > 1 ns in monolayer M 0 S 2 


To compare with the monolayer M 0 S 2 , the polarization resolved PL measurement has 


also been done on the Bernal-stacked (see Fig. 43 i) bilayer M 0 S 2 . Schematic diagram of the 
valleys at the K and K' points with the corresponding spin configurations of the carriers, 
and the diagram for the optical transitions are shown in Figs. |43[e) and (f), respectively. 


Fig. 44 d) and (e) show the observed PL spectra and the measured helicity for the bi¬ 
layer, respectively. The observed low helicity r] exp (bilayer)= 0.25±0.05 for the A exciton 
is in accordance with the hole-spin lifetime ( t v ) of a few hundreds of femtoseconds (from 
the observed 20 times less QY) [2§] , The hole-spin lifetime for the bilayer, therefore, is 
more than three orders of magnitude less compared to the monolayer M 0 S 2 , say by r*4 /. 
The estimated is attributed to the inversion symmetry of the bilayer M 0 S 2 . The <j_ 


polarization couples charge carriers only from one valley K (see Fig. 44 Q, and excites an 


exciton with an electron having spin down and a hole of spin up at the K point. Now the 
intra-valley scattering involves the spin flip of hole which is forbidden by the spin-filtered 
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FIG. 44. (color online)-Left column represents the PL spectra for (a,g and j) single layer and 
(d) for bilayer M 0 S 2 (at 14K). Here the arrows indicate the laser energy used in the experiments. 
The left-handed circularly polarized light (<r_) is used for all the excitations. The middle column 
represent the corresponding observed helicity in PL emissions. The right column represents the 
schematics of the corresponding absorption and emission of lights. Taken from Ref. [ 29] , 


energy gap of ~ 160 meV in absence of magnetic impurity scatterings. The inter-valley 
scattering (from K to K*) involves large momentum transfer with the requirement of the 
spin-flip as well, and hence this channel is also forbidden. Because of the absence of the 
inter as well as intra-valley scattering, the high valley lifetime is observed for monolayer 
M 0 S 2 (29]. Bilayer M 0 S 2 has inversion symmetry which ensures that the Tq(k) = £^(-k). 
The both TRS and the inversion symmetry implies that at each k point the spin degeneracy 
remains irrespective of valley index i.e. E^(k) = £q(k). Hence, for the bilayer the spin and 
the valley index are not coupled and with the cr_ excitation populations of both the valleys 


(K and K') having electron of down-spin and hole of up-spin (see Fig. 44’) are achieved, 
which allows intravalley hole-spin relaxation process via Elliot-Yafet processes ra leading 
to the low observed helicity. 
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FIG. 45. (color online)-PL spectra observed for single layer M 0 S 2 (a) on h-BN substrate and (c) 
on Si02/Si substrate; and the corresponding observed helicity in (b) and (d), respectively. Taken 
from Ref. [ 29] . 


The zero helicity is observed from monolayer M 0 S 2 with the cr_ polarization under the 


laser excitation of 2.33 eV (532 11111 ), as shown in Figs. 44 g) and (h). Since, the probing 
has been done with the higher energy than the resonance energy of A and B excitons, it 


populates the charge carriers at both the valleys, as shown in Fig. 44 i). In this case, no par¬ 
ticular spin-valley selection is occurred and hence the observed zero spin-valley polarization. 


Figs. 44 j) and (k) show the PL spectrum and the helicity of monolayer M 0 S 2 with the laser 
of 2.09 eV energy (594 11111 ), respectively. Because of the resonance with the B excitons, 
the observed helicity is ~ 1, whereas the helicity for the A exciton is zero. I 11 this case, the 
B-exciton relaxes through non-radiative channel to the A-exciton states (phonon-assisted 


process) and populates again both the valleys, as shown in Fig. 44 (1). The polarization 
resolved PL spectra for the monolayer MoS 2 deposited on BN and Si0 2 /Si substrates are 


shown in Figs. |45[a) and (c) and the corresponding measured helicity in Figs. 45 b) and 
(d), respectively. The intensity for BN substrate is ~ 10 times larger than that on SiC^/Si, 
whereas the helicity is almost same for both the substrates. The robustness of the observed 
valley selective CD indicates how good the charge carriers preserve the information of the 
valley index and the spin index with them in their lifetimes i.e. before recombination. 
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FIG. 46. (color online)-Temperature dependent variations of the circular polarizations observed 
in PL spectra for single layer M 0 S 2 . The experimental data have been fitted by considering the 
inter-valley scattering due to the acoustic phonons. Taken from Ref. ]9lj. 


The temperature dependent hclicity of the PL emission spectra for monolayer M 0 S 2 is 


shown in Fig. 46 The circular polarization remains almost constant in the range of 0 to 
90K. With the increasing temperature, the observed hclicity decreases. As the r] exp oc r v , 
with the increase of temperature the phonon populations and hence the scattering via them 
also increases by providing the required large momentum transfer from K to K' point in BZ. 
The inter-valley scattering rate is t^ 1 oc e~ Eph -/ KBT , where E p h is the phonon energy at the 
K-point. The extracted value of the phonon frequency by fitting the observed PL hclicity 
is given by E p h & 240 cm -1 in accordance with the acoustic phonon near the K-point of 
the monolayer M 0 S 2 [63]. Although, the fitting is in agreement with the observed hclicity 
above 90K, more experimental as well as theoretical studies are needed to explore the full 
understanding of the underlying physics. 
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9.4 . Gate-voltage dependent absorption and photoluminescence of the A exciton 


complexes: observation of tightly bound trions 


A few absorption spectrum for the monolayer M 0 S 2 at different back gate voltages in the 
range of -100 to +80 V using 532 nm laser line are shown in Fig. [47^ a). As we discussed, the 
monolayer M 0 S 2 device shows n-type transfer characteristics at ambient conditions because 
of the Fermi level pining. Therefore, the device operated at -100V is having almost the 
undoped channel. Here we discuss the absorption and PL spectra for the AEC only. The 
absorption spectrum at -80V is dominated by an intense neutral A-exciton and an little 
shoulder due to the negatively charged A-exciton. This charged exciton consisting of one 
hole of up-spin and two electrons of opposite spins is known as trion (A - ) [30]. As the 
negative gate voltage decreases, the A~ evolves and around +70V, the spectrum is domi¬ 
nated completely by A~ exciton. The neutral A-exciton peak gradually diminishes with the 
decrease of the negative gate voltages, and after a certain voltage (0V) it is hard to track the 
peak. The peak positions of these two peaks are plotted versus back gate voltages (Fermi 
energy in lower axis) in Fig. |4T](b) , and their difference (c +4 — to a- ) versus Fermi energy (gate 
voltages in upper axis) is plotted in Fig. |47](c). The overall reduction of the intensity for 
the absorption spectra is due to the Pauli blocking associated with the increasing electrons 
in the conduction band |30j . In the right panel of Fig. |47](a), the PL spectra at different 
gate voltages is shown using 532 nm laser line. The neutral A-exciton behaves similarly 
as in absorption spectra i.e. after a positive gate voltage it is turned off. The feature of 
the blueshifts of the neutral A-exciton and the almost constant feature of the A~-exciton 
with gate voltages are attributed to the combined effect of the many-body interactions and 
the Pauli blocking [30]. The neutral A-exciton can be thought of an ionized trion such 
that the splitting between the peaks is given by loa — <+ 4 - = Ea- + Ep\ where, Ea- is the 
binding energy of the trion A~. At Ep= 0, Ea- is the separation between the two peaks 
and the extracted value by fitting the linear relation (see Fig. 47 2) is ~ 18.0+1.5 meV. The 


separation uja — oja- behaves linearly with the Fermi energy ( Ep ), and is in accordance with 
the earlier observations of trions for the quantum wells [ 95] . 


Figs. |48)(a) and (b) show the room temperature PL spectra (with broadened features) 
and the corresponding gate voltage dependent measured PL intensity, respectively. The PL 
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FIG. 47. (color online)-(a) Back-gate voltage dependent absorption spectra (left panel) and the 
PL spectra (right panel) measured at 10K. The power-law has been used to fit (blue solid lines) 
the data for the absorption spectra and Lorentzian fits has been done for PL spectra. Here the 
dashed blue lines are guide to the eye. The arrows in PL spectra indicate the shift of the trion 
with the gating, (b) The extracted (from absorption spectra) peak positions of trion (oja-) and 
the neutral exciton (oja) as a function of gate voltages (upper scale) and the corresponding Fermi 
energy (bottom scale); and the difference u>a — hi ( c ) • The E p dependence linear fit has been 
done to the data. The inset in (c) shows the schematics of the breaking of the trion into a neutral 
exciton and an electron at the Fermi level. Here the offset gate voltage =-107+6 V for E^=0. 
Taken from Ref. m- 


intensity for the neutral A-exciton decreases with gating by almost two orders of magni¬ 
tude, whereas for the trion it remains constant in the range of -100 to +80 V. The decreased 
PL intensity of the neutral A-exciton is attributed to the spectral weight reduction and the 


transfer of it to the trion with the doping as for the absorption spectra at 10K (see Fig. 47 1 ). 
This high tunability of the PL intensity, which is not observed at 10K, is attributed to the 
thermally activated carrier populations from the trion state. 
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FIG. 48. (color online)-(a) The gating voltage dependent PL spectra measured at 300K. (b) The 
integrated intensity of both the neutral exciton and the trion and their total intensity (left scale). 
Right scale represents the transfer characteristics. Taken from Ref. [30]. 


The gating voltage dependent optical response might be associated with the metal- 
insulator transition (MIT) as observed for the two-dimensional electron gas (2DEGs) [9B]. 
At ~ -100V, the channel is undoped (T=10K) and hence the optically excited excitons 
dominate the absorption spectra. Now as we increase the doping, for positive gate voltages 
the conduction band is populated with the free carriers and the neutral A-exciton combines 
with them to form trion. Hence, at higher positive gate voltages, the absorption spectra is 


dominated completely by the trions (see Fig. 47 r). As we depletes the channel more and 
more, the decrease in the number of carriers leads to the much less screening of the long range 
Coulomb interactions due to the presence of the extrinsic disorders (ionized) in the sample. 
This implies an increase in potential fluctuations in the sample which in turn localize the 
electrons. This feedback mechanisms results in non-linear screening and hence, can reduce 
the MIT region to a smaller gate voltage. Recently, it has been shown that the monolayer 
M 0 S 2 indeed goes to MIT at the top gate voltage of ~ 2.2V (the carrier concentration ~ lx 


10 13 cm -2 ) below 80K [31 j, as shown in Fig. 49 If 2DEG was uniformly distributed instead 
of localized, we should not observe the improvement of the spectral weight of the neutral 
A-exciton peak in absorption spectra. Since localized electrons do not play an effective 
role to the screening, the more localization means more prominent neutral A-exciton peak. 


67 








Now, if the charge carriers would form the patches of inhomogeneous carrier density, the 
associated tuning of the PL intensity should not show any on-set. In contrast, experimental 
results for the PL intensity show an on-state (threshold) with the gate voltage (see Fig. [48 Jd) , 
which implies the possibility of the localization of electrons (Wigner crystallization) instead 
of forming small patches [96]. If the localized electrons were effective in the screening, then 
we should not observed the enhancement of the A-exciton in absorption spectra; and hence, 
we can say that this MIT may be due to the Mott-Hubbard type instead of Anderson-type. 
More experimental works are required to explore the physics in the insulating region and 
the corresponding theoretical supports to explain the observed MIT in the monolayer M 0 S 2 . 




FIG. 49. (color online)-(a) The top gate (Vt g ) dependent conductance (G) of single layer M 0 S 2 for 
various temperatures. The G shows the thermally activated behavior for low Vt g and it decreases 
with decreasing temperature. Single layer M 0 S 2 shows a transition from an insulating regime to 
a metallic regime above Vtg of ~ 1-2 for the decreased temperatures such that the G increases 
with decreasing temperature and this cross over is clear for the lower temperature curves, (b) The 
temperature dependent G of the devices for different values of the electron density ( 112 / 3 )- Taken 
from Ref. El- 


10. CONCLUSIONS 

In conclusion, single and a few layer M 0 S 2 devices show reasonably high carrier mobility 
(~ 700 cm 2 /V-sec) and on-off ratio ~ 10 8 . The increment of the high temperature (100-300 
K) mobility has been engineered by quenching the homopolar and other phonon scattering 



















processes with the help of dielectric coating. The observed robustness of the valley selective 
circular dichrosim suggest that we should explore experimentally valley Hall effect and mag¬ 
netism in single-layer M 0 S 2 without applying magnetic field. The observation of negatively 
charged excitons or trions and the associated metal-insulator-transition at low temperature 
can be an indication of Wigner-crystal formation in 2-D M 0 S 2 material, which needs to be 
further explored. Along with that, the measured dimensionless parameter (ratio of Coulomb 
potential energy to kinetic energy) r s « 60 make the M 0 S 2 as an ideal system for further 
exploration of many-body phenomena. In addition, other dichalcogenide materials like WS 2 , 
MoSe 2 and WSe 2 [97j are exciting to explore the novel physics in coming years. 
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